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FOREWORD 


The  work  described  in  this  report  is  based  on  an  invited  lecture  series  by  the  Advisory 
Group  for  Aerospace  Research  and  Development  (AGARD)  to  be  given  in  June  1994  at  the  Von 
Karman  Institute  (VKI)  in  Belgium  and  the  Middle  East  Technical  University  in  Ankara,  Turkey. 
While  the  lecture  series  is  generic  in  nature,  it  focuses  on  a  Class  of  Aeroprediction  codes  that 
are  semiempirical  in  nature  (a  mixture  of  analytical  and  empirical  methods).  By  design  of  the 
invitation,  the  report  also  emphasizes  many  of  the  new  methods  the  author  has  developed  over 
his  thirty  year  career  at  NSWCDD.  These  lectufes  are  being  documented  in  the  form  of  a 
technical  report  in  order  that  those  who  are  not  able  to  attend  the  presentations  will  have  easier 
access  to  the  material 

While  the  new  methods  are  the  authors  principally,  implementation  of  these  methods  was 
done  primarily  by  coworkers.  As  a  result,  a  lengthy  acknowledgement  is  given  at  the  end  of 
the  report  which  includes  not  only  the  technical  personnel  but  sponsors  as  well.  This  particular 
report  was  supported  by  the  Office  of  Naval  Research  (Dave  Siegel)  and  more  specifically,  the 
Surface  Launched  Weapons  Technology  Program  managed  at  NSWCDD  by  Robin  Staton  and 
the  Air  Launched  Weapons  Technology  Program  managed  at  NAWC/China  Lake  by  Tom 
Loftus.  Appreciation  is  expressed  to  these  individuals  for  their  support. 

Approved  by: 

DAVID  S.  MALYEVAC,  Deputy  Department  Head 

Weapons  System  Department 
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ABSTRACT 


This  report  discusses  the  pros  and  cons  of  numerical,  semiempirical  and  empirical 
aeroprcdiction  codes  and  lists  many  state-of-the-art  codes  in  use  today.  It  then  summarizes  many 
of  the  more  popular  approximate  analytical  methods  used  in  State-of-the-Art  (SOTA) 
semiempirical  aeroprediction  codes.  It  also  summarizes  some  recent  new  nonlinear 
semiempirical  methods  that  allow  more  accurate  calculation  of  static  aerodynamics  on  complete 
missile  configurations  to  higher  angles  of  attack.  Results  of  static  aerodynamic  calculations  on 
complete  missile  configurations  compared  to  wind  tunnel  data  are  shown  for  several 
configurations  at  various  flight  conditions.  Calculations  show  the  new  nonlinear  methods  being 
far  superior  to  some  of  the  former  linear  technology  when  used  at  angles  of  attack  greater  than 
about  15  degrees. 
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1.0  INTRODUCTION  AND  BACKGROUND 


1 . 1  USES  FOR  AERODYNAMICS 


Aerodynamics  are  required  throughout  the  design  process  of  any  flight  vehicle.  These 
aerodynamics  are  used  for  flight  performance  estimates  including  range,  maneuverability, 
miss  distance,  and  stability  analysis.  In  addition,  they  are  used  for  structural  analysis 
including  material  requirements  and  selection,  structural  member  thicknesses  required  to 
withstand  the  loads,  and  as  inputs  for  heat  transfer  or  ablation  analysis  (Table  1-1). 

Generally,  an  interactive  design  process  occurs  between  the  aerodynamicist,  the  structural 
designer,  and  the  flight  dynamicist  to  arrive  at  a  configuration  that  meets  some  set  of  desired 
launcher  constraints  and  performance  requirements  given  a  warhead  and  possibly  a  guidance 
system  as  well. 


TABLE  1-1.  WHAT  AERODYNAMICS  ARE  USED  FOR 


Flight  Dynamics 

Structures 

o  Range  Computation 

o  Engagement  of  Target  and  Miss  Distance 

o  Maneuverability  Estimates 

o  Any  Trajectory  Analysis  (3  DOF,  5  DOF, 

6  DOF)* 

o  Loads  (Pressure) 

o  Aeroheating  (Inputs  to 

Heat  Transfer  Codes) 
o  Ablation  Analysis  Inputs 

DOF  =  Degree  of  Freedom 


Prior  to  1971,  the  tactical  weapons  aerodynamicist  could  do  one  of  three  things  to 
obtain  aerodynamics.  The  individual  could  perform  flight  tests  of  a  full-scale  configuration; 
or  design,  build,  and  test  a  wind  tunnel  model  over  the  flight  range  of  interest;  or  finally, 
utilize  existing  handbooks,  wind  tunnel  data  reports,  and  theoretical  analysis  to  estimate 
empirically  the  aerodynamics  of  a  given  configuration. 

The  first  two  approaches  were  often  more  costly,  time  consuming,  and  accurate  than 
needed  in  the  preliminary  design  stages,  whereas  the  latter  approach  was  more  time 
consuming  than  desired  but  also  had  no  general  accuracy  assessment. 

A  fourth  alternative  (which  did  not  exist  prior  to  1971),  to  compute  aerodynamics  on 
a  complete  configuration  over  the  Mach  number  and  angle  of  attack  range  of  interest,  is  to 
have  a  general  computer  program  to  perform  such  a  task.  There  are  three  alternative 
theoretical  approaches  to  develop  such  a  code  (see  Table  1-2).  The  first  of  these  is  solution 
of  the  full  Navier  Stokes  equations.  The  only  assumptions  associated  with  this  set  of 
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equations  is  continuum  flow  (that  is  the  flowfield  region  is  not  sparsely  populated  with  air 
molecules  such  as  at  altitudes  greater  than  about  200  to  250  thousand  ft)  and  the  turbulence 
model  selected.  A  second  theoretical  alternative  is  to  assume  the  viscous  flow  region  lies  in 
a  thin  layer  near  the  body  and  thus  solution  of  the  Navier  Stokes  equations  can  be  reduced  to 
that  of  an  inviscid  flowfield  plus  a  thin  boundary  layer  near  the  surface.  This,  combined 
with  empirical  estimates  of  base  drag  and  other  protuberance  aerodynamics,  gives  a  complete 
set  of  aerodynamics  for  the  configuration  of  interest.  A  third  theoretical  alternative  is  to 
assume  the  body  perturbs  the  flowfield  only  slightly  and  then  to  make  appropriate 
approximations  to  the  Euler  and  Boundary  Layer  Equations.  These  approximate  theories  are 
then  combined  with  other  theoretical  approaches  and  empirical  data  for  the  complete 
aerodynamics  code. 


TABLE  1-2.  HOW  WE  GET  AERODYNAMICS 


1. 

Wind  Tunnel,  Free  Flight  Data,  Ballistic  Range 

2. 

Empirical  Estimates:  Wind  Tunnel  Reports,  Handbooks,  Experience,  etc. 

3. 

Aeroprediction  Codes 

A. 

Navier  Stokes  —  Continuum  Flow 

B. 

Euler  Equations  +  Boundary  Layer  --  inviscid  outer  layer  +  thin 
viscous  layer  near  surface  +  some  empirical  techniques 

C. 

Approximations  to  Euler  and  Boundary  Layer  Equations  + 

Empirical  Techniques 

There  are  several  uses  that  can  drive  the  type  of  theory  chosen  for  the  aeroprediction 
code.  These  are  listed  in  Table  1-3.  For  example,  if  missile  synthesis  is  being  performed 
where  a  very  large  number  of  configurations  are  investigated  to  conduct  top  level  trade 
studies  involving  engine  types,  warhead  types,  material  requirements,  etc.  as  a  function  of 
range,  maneuverability,  or  response  time,  then  it  is  desirable  to  have  an  easy  to  use,  robust, 
and  computationally  fast  code.  At  the  same  time,  accuracy  may  be  sacrificed  to  achieve 
these  goals. 

After  a  missile  synthesis  of  a  large  number  of  concepts  has  been  conducted,  generally 
several  of  these  concepts  are  taken  a  step  further  in  the  design  process.  Here,  structural 
layouts,  packaging  of  all  components,  and  better  definition  of  weights  are  typical 
requirements  that  allow  improved  estimates  of  range,  maneuverability,  and  preliminary  miss 
distance.  This  means  that  the  aerodynamic  code  requirements  need  a  blend  of  robustness, 
ease  of  use,  and  accuiacy  while  still  being  computationally  cost  effective.  Accuracies  in 
aerodynamics  of  10  percent  or  so  are  generally  expected. 

Finally,  one  or  two  configurations  are  selected  for  more  detailed  performance 
estimates.  This  means  accuracy  in  the  aerodynamics  estimates  of  better  than  5  percent  in 
most  cases.  Each  of  the  three  design  levels  discussed  require  different  levels  of  accuracy, 
computational  speed,  and  robustness  and,  therefore,  aid  in  the  choice  of  the  level  of 
theoretical  complexity  needed  to  meet  the  requirements. 
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TABLE  1-3.  AERODYNAMIC  CODE  REQUIREMENTS  AND  USES  IN  VARIOUS 

MISSILE  DESIGN  STAGES 


Design  Stage 

Aero  Code  Design 
Requirements 

Trade  Studies 
(Typical) 

Aerodynamics 

Uses 

Missile  Synthesis 

Robustness 

Ease  to  Use 

Minimal  Input 
Parameters 

Extremely  Fast 
Computationally 

25  percent  Accuracy 

Engine  Types 

Warhead  Types 
Material 

Requirements 

Typical  Weights 
Guidance  Types 
Airframe  Control 

Type 

Range 

Maneuverability 
Response  Time 

Missile  Preliminary 
Design 

Blend  of  Robustness, 

Ease  of  Use,  and 
Accuracy 

Fast  Computationally 

10  percent  Accuracy 

Structural  Layout 
(Material, 

Thickness,  etc.) 

Aero  Shape  vs. 
Engineering  and 
Guidance  Size 

Hot  vs.  Cold 

Structure 

Range 

Maneuverability 
Miss  Distance 
(3  DOF) 

Structural  Design 

Detailed  Design 
and  Problem 

Solving  (or 

Analysis  Codes) 

Accuracy  (<5  percent) 
Computationally 
Affordable 

User  Friendliness  and 
Robustness  Still 
Important 

Detailed  Structural 
Design  Including 
Material  Selection 
Investigating 

Critical  Problem 
Areas 

Range 

Maneuverability 
Miss  Distance 
(6  DOF) 

Structural  Design 

To  meet  the  theoretical  aerodynamics  computer  code  needs,  the  Navy  began 
developing  such  a  code  in  1971  based  on  the  3C  approach  of  Table  1-2.  This  code  falls  into 
the  second  category  of  Table  1-3.  Since  the  first  version  of  the  NSWCDD  Aeroprediction 
code  was  released,  there  have  been  four  versions  produced  since  that  time. 

Each  of  these  versions  attempted  to  meet  the  requirements  as  seen  by  the  tactical 
weapons  community.  The  first  version  was  for  general-shaped  bodies  alone.1  It  was  the  first 
such  weapons  code  known  that  combined  a  good  ir.  x  of  accuracy  in  aerodynamic 
computations,  ease  of  use  and  computational  time.  It  is  believed  that  this  mix  led  to  the 
code’s  initial  popularity  and  requests  for  additional  capability.  In  1974, 2>  3  the  code  was 
extended  to  allow  up  to  two  sets  of  lifting  surfaces  in  the  computational  process.  In  1977 ,4, 5 
dynamic  aerodynamic  derivatives  were  added  to  the  code’s  capability.  In  1981,  the  code 
extended  the  Mach  number  range  up  to  eight  and  added  high  angle-of-attack  capability  for  a 
narrow  range  of  configurations6’7.  Finally,  the  last  version  of  the  code  extended  the  Mach 
number  range  higher  to  include  real  gas  effects,  added  new  nonlinear  lift  methodology  for 
wings  and  interference  effects,  and  developed  an  improved  base  drag  methodology.*’9 
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This  report  will  serve  several  purposes.  First,  a  review  of  the  state-of-the-art  (SOT A) 
semiempirical  aerodynamic  prediction  codes  will  be  given.  Second,  a  review  of  some  of  the 
more  useful  approximate  theoretical  methods  will  be  made.  These  methods  are  conventional 
and  have  been  in  use  for  many  years.  Third,  a  more  detailed  review  of  the  new  nonlinear 
aerodynamic  methods  introduced  over  the  past  3  years  into  the  fifth  version  of  the 
Aeroprediction  Code  (AP93)  will  be  given.  Finally,  a  comparison  of  the  static  aerodynamics 
using  experiment,  the  AP93  and  the  older  version  of  the  Aeroprediction  Code  (AP81)  will  be 
made  on  several  complete  missile  configurations. 


1.2  TYPES  OF  AEROPREDICTION  CODES 


Aeroprediction  Codes  will  be  defined  and  broken  down  into  three  classes.  These 
classes  are  empirical,  semiempirical,  and  numerical  codes.  The  empirical  codes  are 
analogous  to  the  codes  used  in  Missile  Synthesis  in  Table  1-3.  The  semiempirical  and  some 
numerical  codes  are  used  primarily  in  the  missile  preliminary  design  stage  of  Table  1-3. 
Finally,  the  numerical  codes  are  the  only  ones  with  the  accuracy  and  capability  to  do  the 
detailed  design  application  as  shown  in  Table  1-3. 

In  terms  of  a  definition,  empirical  codes  typically  calculate  aerodynamics  by  a  series 
of  simple  formulas  that  have  been  approximated  based  on  data  fits.  Typically,  these  codes 
can  be  implemented  on  a  hand  calculator  in  many  cases  and  are  the  most  simplistic  and  least 
accurate  of  the  code  classes. 

The  semiempirical  codes  typically  attempt  to  calculate  a  force  or  moment  using 
approximations  to  the  exact  equations  of  motion.  When  this  approach  fails  (such  as  at  higher 
angles  of  attack),  empirical  estimates  or  methods  are  used.  This  blend  of  approximate 
theories  and  empirical  estimates  is  why  this  class  of  codes  is  termed  semiempirical.  The 
semiempirical  codes,  in  contrast  to  the  empirical  codes,  generally  will  calculate  pressure 
distribution  on  the  body  and  lifting  surfaces.  It  is  this  blend  of  theory  with  the  empirical 
estimates  that  allows  the  semiempirical  codes  to  improve  accuracy  over  the  empirical  codes. 

The  third  class  of  codes  is  called  numerical.  These  codes  will  define  a  grid  around 
the  configuration  that  is  composed  of  points  in  two  or  three  dimensions.  Numerical 
techniques  are  then  employed  to  solve  the  equations  of  motion  at  all  grid  points  in  the  flow 
field  that  is  bounded  by  the  body  and  shock  or  body  and  outer  boundary  of  the  flow  if  the 
Mach  number  is  subsonic.  Numerical  Codes  are  generally  based  on  the  linearized  or  full 
potential  equations  of  motion,  the  full  Euler  equations  or  the  full  or  reduced  level  of  Navier 
Stokes  equations.  If  the  potential  or  Euler  equations  are  used,  other  methods  (such  as 
boundary  layer  equations)  must  be  used  for  skin  friction.  Also,  empirical  estimates  are  used 
for  base  drag.  Hence,  even  though  these  codes  are  numerical,  in  most  cases  to  get  complete 
forces  and  moments  on  a  configuration,  the  use  of  some  empirical  data  will  be  necessary. 
Also,  if  the  potential  equations  are  solved  in  a  numerical  form,  the  accuracy  is  similar  to  the 
semiempirical  codes.  The  only  difference  between  the  two  is  that  the  semiempirical  codes 
seek  pressure  distributions  on  the  body  and  wings  without  solving  the  entire  flowfield.  This 
saves  a  tremendous  amount  of  computational  time. 
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A  final  point  worthy  of  discussion  are  the  assumptions  inherent  in  each  level  of 
theory.  These  assumptions  are  given  as  a  function  of  the  theoretical  approach  in  Table  1-4. 
Upon  examination  of  Table  1-4,  the  level  of  code  sophistication,  computational  time,  overall 
cost  and  accuracy  goes  down  in  going  from  the  top  to  the  bottom  of  the  table. 

One  way  to  try  to  compare  the  level  of  sophistication  versus  accuracy,  and  the  cost  of 
the  various  codes,  is  through  the  examination  of  die  total  cost  to  obtain  a  set  of 
aerodynamics.  To  do  this,  Table  1-5,  which  compares  the  educational,  computer,  and 
computational  time  requirements  of  the  various  Aeroprediction  Codes  in  use  at  NSWCDD 
has  been  prepared.  Referring  to  Table  1-5,  the  level  of  sophistication  increases  in  going 
from  top  to  bottom  of  the  Table.  For  example,  the  MAIR  Code  is  close  to  an  empirical  code 
but  it  does  have  some  theory  included  so  that  it  would  be  in  the  class  of  semiempirical  codes. 
The  Missile  III,  Aeroprediction  versions  81  and  93,  HABP,  and  missile  DATCOM,  are  all 
semiempirical  codes.  NANC  and  BODHEAT  are  primarily  numerical  codes  based  on 
approximations  to  the  Euler  and  Boundary  Layer  equations.  SWINT/ZEUS,  CFL3DE  and 
GASP,  of  course,  are  all  numerical  codes.  The  Aeroprediction  81/93,  SWINT/ZEUS, 

MAIR,  NANC,  and  BODHEAT  were  all  developed  at  NSWCDD.  The  Missile  III  was 
developed  by  Nielsen  Engineering  and  Research  (NEAR),  HABP  and  Missile  DATCOM  by 
McDonnel  Douglas  of  St.  Louis,  and  the  Navier  Stokes  Codes  were  developed  jointly  by 
NASA/LRC  and  VPI. 

Included  in  Table  1-5  is  the  time  required  to  leam  how  to  use  the  code,  the  set-up 
time  for  a  typical  geometry,  and  the  computer  time  for  the  one  case  referenced  to  the  same 
computer  (CDC  865).  Also  shown  are  other  criteria  including  typical  educational  level  of 
the  user  as  well  as  the  size  of  the  computer  required.  To  get  the  total  cost  of  using  a  code,  it 
is  necessary  to  add  the  manpower  set-up  time  to  the  computer  cost  and  prorate  the  training 
time  over  some  nominal  expected  usage.  Experience  has  shown  that  most  project  and 
program  managers  are  willing  to  pay  the  costs  of  SWINT/ZEUS  type  codes  and  any  above 
that  in  Table  1-5.  However,  the  cost  and  requirements  of  the  full  Navier  Stokes  codes  must 
come  down  substantially  before  they  will  be  used  on  a  routine  basis  for  design.  This  means 
much  additional  research  as  well  as  advancements  in  computer  speed  are  still  needed  in  this 
area. 

To  illustrate  this  point,  a  particular  example  was  chosen  for  cost  comparisons.  The 
example  is  to  develop  a  set  of  trim  aerodynamics  on  a  typical  missile  configuration  to  be 
used  as  an  input  to  a  three-degree-of-freedom  (3  DOF)  flight  simulation  model.  This 
example  is  quite  typical  of  what  an  empirical  or  semiempirical  code  would  be  used  for.  By 
definition,  trim  is  that  combination  of  angles  of  attack  (a’s)  and  control  deflections  (5’s)  that 
give  zero  pitching  moment  about  the  vehicle  center  of  gravity.  To  determine  the  (a,  5)  map 
as  a  function  of  Mach  number,  one  must  compute  the  static  aerodynamics  over  enough  a,  5, 
M  conditions  so  the  flight  envelope  will  be  covered.  Also,  it  will  be  assumed  that  the 
missile  is  a  surface  launched,  tail  control,  cruciform  fin  configuration  which  has  a  Mach 
range  of  0  to  4,  angle  of  attack  range  of  0  to  30°,  control  deflection  of  0  to  20°,  and  attitude 
0  to  80,000  feet.  These  conditions  are  reasonable  for  many  of  the  worlds  missiles.  To  cover 
the  flight  envelope,  7  Mach  numbers,  5  a’s  and  5  5’s  are  assumed.  This  gives  a  total  of 
7x5x5  =  175  cases.  Furthermore,  skin  friction  varies  with  attitude  so  5  altitudes  will  be 
chosen,  giving  a  total  of  180  cases  for  which  aerodynamics  are  to  be  computed  on  a  single 
configuration. 
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TABLE  1-4.  ASSUMPTIONS  OF  FLOW  FIELD  EQUATIONS 

1.  Full  Navier  Stokes  (high  angle  of  attack) 

A.  Continuum  Flow 

B.  Turbulence  Model 

2.  Thin  Layer  Navier  Stokes  (moderate  separation) 

A.  Neglect  Streamwise  and  Circumferential  Gradients  of  Stress  Terms 

B.  Turbulence  Model 

C.  Continuum  Flow 

3.  Parabolized  Navier  Stokes  (small  separation) 

A.  Steady  State 

B.  Neglects  Streamwise  Viscous  Gradient 

C.  Approximate  Streamwise  Pressure  Gradient  in  Subsonic  Portion  of  Flow  Near  Surface 

D.  Turbulence  Model 

E.  Continuum  Flow 

4.  Euler  Equations  +  Boundary  Layer  (small  separation) 

A.  Viscous  Region  Confined  to  Thin  Region  Near  Body  Surface 

B.  Large  Reynold’s  Number 

C.  Neglect  Streamwise  Gradients  of  Stress  Terms 

D.  Neglect  Normal  Pressure  Gradient 

E.  Turbulence  Model 

F.  Continuum  Flow 

5.  Euler  Equations 

A.  Neglect  all  Viscous  Terms 

B.  Continuum  Flow 

6.  Full  Potential  Equations 

A.  Neglect  all  Viscous  Terms 

B.  Flow  is  Isentropic  (no  shock  waves) 

C.  Continuum  Flow 

7.  Linearized  Potential  Equations 

A.  Neglect  all  Viscous  Terms 

B.  Flow  is  Isentropic  (no  shock  waves) 

C.  Body  Creates  Small  Disturbances  in  Flowfield 

D.  Continuum  Flow 

8.  Theoretical  Approximations 

A.  Certain  Other  Simplifications  to  Euler,  Potential  Equations,  or  Boundary  Layer  Equations 

B.  Continuum  Flow 

9.  Empirical  Data  Base 

A.  Data  Base  Covers  Vehicles  and  Flight  Regime  of  Interest 
_ B. _ Enough  Data  is  Available  to  do  Good  Interpolations _ 
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TABLE  1-5.  EDUCATIONAL  AND  TIME  REQUIREMENTS  FOR  AEROPREDICTION  CODES  IN  USE 


AT  NSWCD] 

D 

Code 

Typical  User 
Educational  Level 

Typical 

Time 

Required  to 
Learn  to 

Use  Code 

Set-Up 

Time 

Computational 
Time  for  1 
Case  (Same 
Computer) 

Computer 

Required 

1. 

MAIR 

Coop,  B.S.,  M.  S., 

Ph.  D 

<  1  wk 

<  1  day 

<  1  second 

n 

2. 

Missile  in 

Coop,  B.S.,  M.  S., 
Ph.  D 

*  1  wlc 

<  1  day 

<  1  second 

n 

3. 

Aeroprediction 

81  and  93 

Coop,  B.S.,  M.  S., 
Ph.  D 

*  1  wk 

<  1  day 

<  1  second 

p.  c. 

■ 

HABP 

B.  S.,  M.  S., 

Ph.  D. 

•  2wk 

<  1  wk 

<  1  second 

Micro 

Vax 

5. 

Missile 

DATCOM 

B.  S.,  M.  S„ 

Ph.  D. 

*  2  wk 

<  1  wk 

<  1  second 

Micro 

Vax 

6. 

NANC 

M.  S.,  Ph.  D. 

m  3  wk 

<  2  wks 

10  seconds 

Vax  CDC 
Super 
Mini 

1 

BODHEAT 

M.  S.,  Ph.  D. 

*  3  wk 

<  1  wk 

10  seconds 

Vax  CDC 
Super 
Mini 

8. 

SWINT/ZEUS 

M.  S.,  Ph.  D. 

*>  1  month 

<  1  month 

1-3  minutes 

Vax  CDC 
Super 
Mini 

9. 

N.S.  (CFL3DE, 
GASP) 

Ph.  D.,  some  M.  S. 

»  months- 

yrs 

m  months 

m  hrs-days 

Cray  or 
Super 
Mini 

Before  costs  of  each  computer  code  can  be  made  for  this  particular  example,  some 
assumptions  must  be  made.  These  assumptions  are  given  in  Table  1-6.  These  assumptions 
are  based  on  NSWCDD  experience  in  using  the  various  aeroprediction  codes.  The  cost  to 
perform  the  set  of  trim  aerodynamics  calculations  using  these  codes  is  shown  in  Figure  1-1. 

It  should  be  noted  that  the  cost  assumes  that  Parabolized  Navier  Stokes  and  Euler  plus 
boundary  layer  are  used  at  subsonic  axial  Mach  number  conditions  although  the  codes  in  use 
at  NSWCDD  are  steady  hyperbolic  marching  solutions  and  will  not  function  where  the  axial 
Mach  number  decreases  to  one.  To  go  to  unsteady  computation  would  require  costs  to  be 
multiplied  by  a  factor  of  at  least  10.  Hence,  the  PNS  and  Euler  plus  B.L.  costs  are  based  on 
steady  flow  of  supersonic  Mach  numbers.  For  a  combination  of  steady  and  unsteady 
computations,  the  cost  of  these  codes  would  probably  be  about  five  times  greater  than  those 
shown  in  Figure  1-1. 
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TABLE  1-6.  ASSUMPTIONS  IN  COST  ESTIMATES  TO  COMPUTE  SET  OF  TRIM 
AERODYNAMICS  WITH  VARIOUS  AEROPREDICTION  CODES 


Estimated  Cost  o 
o 
o 

o 

o 


Cray  n  Computer  at  $500/HR 
Engineer  Time  =  1  1  OK/work  year 

Engineer  is  assumed  to  know  how  to  use  codes  so  no  training 
time  is  involved. 

Need  enough  resolution  in  grid  size  to  predict  skin  friction  drag 
Wind  Tunnel  (W/T)  includes  model  and  test  cost 


CODE 

FNS 

TLNS 

PNS 

EULER  +  BL  +  B.D. 
AEROPREDICTION 


SET  UP  TIME 
5  Weeks 
5  Weeks 
2-5  Weeks 
2  Weeks 
0.5  Day 


COMPUTER  TIME 
20  Hours 
17  Hours 
12  Minutes 
1.5  Minutes 
1.0  Seconds 


There  are  several  points  worthy  of  note  in  analyzing  Figure  1-1.  First,  for  practical 
routine  computations,  Full  Navier  Stokes  and  Thin  Layer  Navier  Stokes  are  beyond  the  cost 
most  program  managers  are  willing  to  pay.  Secondly,  they  are  even  beyond  the  wind  tunnel 
cost  to  obtain  comparable  aerodynamics.  Thirdly,  steady  PNS,  steady  Euler  plus  boundary 
layer,  and  semiempirical  (Aeroprediction)  are  all  within  most  allowable  aerodynamics 
budgets.  Going  to  unsteady  computations  for  subsonic  axial  Mach  numbers  makes  the  cost 
requirements  much  higher  and  may  not  be  affordable  and  robust  to  cover  the  entire  flight 
regime. 

A  second  way  of  comparing  aerodynamic  computations  is  the  total  time  it  takes  to  get 
the  complete  set  of  computations  performed.  These  results  are  estimated,  again  based  on 
NSWCDD  experience,  and  shown  in  Figure  1-2.  Again,  the  same  caveat,  with  respect  to  the 
PNS  and  Euler  Codes,  applies  here  as  to  Figure  1-1.  For  most  development  programs,  the 
semiempirical  codes  obviously  have  the  most  desirable  tum-around-dme  (TAT).  The  Euler 
and  PNS  are  marginal  and  experimental  and  Navier-Stokes  (N-S)  and  Thin  Layer  Navier- 
Stokes  (TLNS)  generally  unacceptable  except  as  long  lead  items.  The  combination  of  cost, 
accuracy,  and  complexity  of  the  various  means  of  computing  aerodynamics  has  led  most 
agencies  to  a  mix  of  the  various  approaches.  The  most  used  codes  still  remain  the 
semiempirical  codes  with  Euler  plus  Boundary  Layer  becoming  more  and  more  prevalent  as 
the  robustness  and  ease  of  use  improves.  Navier  Stokes  and  Thin  Layer  Navier  Stokes  are 
used  for  specialized  problems  or  a  few  validation  cases  of  other  codes;  much  work  is  still 
needed  to  improve  user  friendliness  for  this  class  of  codes.  Wind  tunnel  data  still  remains 
the  most  reliable  but  time  consuming  method  to  obtain  Aerodynamics. 
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CODE  COMPLEXITY 


FIGURE  1-1.  ESTIMATED  COST  TO  OBTAIN  SET  OF  TRIM  AERODYNAMICS 
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TIME 

(MONTHS) 


FIGURE  1-2.  ESTIMATED  TIME  TO  GENERATE  SET  OF  TRIM  AERODYNAMICS 
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Codes  in  Use 

Tables  1-7  through  1-12  list  many  of  the  aerodynamic  codes  in  use  today.  These  tables 
were  taken  primarily  from  Lacau,10  with  some  of  the  more  important  codes  developed  since 
1988  added  to  the  tables.  References  11  -  57  are  the  references  one  could  use  for  each  of 
the  codes  listed  in  Tables  1-7  through  1-12. 10  (Note  that  some  of  the  references  are 
incomplete  but  were  taken  directly  from  Reference  11.)  The  tables  are  given  in  a  similar 
fashion  to  the  theoretical  assumptions  given  in  Table  1-4.  All  the  previous  comments,  with 
respect  to  the  general  types  of  codes,  should  also  be  kept  in  mind  when  considering  a 
particular  type  of  code  to  accomplish  a  given  task.  No  attempt  will  be  made  here  to  expose 
the  good  and  bad  points  of  each  of  the  codes.  This  would  require  a  personal  knowledge  in 
terms  of  usage  on  a  set  of  configurations  for  a  given  problem(s).  This  obviously  would  be 
prohibitively  expensive  and  time  consuming.  The  best  that  one  can  generally  hope  for,  is  a 
comparison  of  a  few  codes  for  a  limited  class  of  conditions. 

This  completes  the  discussion  on  the  state-of-the-art  in  aerodynamic  codes  and  the 
various  means  to  obtain  aerodynamics.  The  bulk  of  the  remainder  of  this  paper  will  be 
directed  at  the  semiempirical  code  known  as  NSWC  Aeroprediction  or  NSWC-AP  as  given 
in  Table  1-7.  To  that  extent,  the  next  section  will  briefly  cover  many  of  the  more  popular 
approximate  theoretical  techniques  used  by  many  of  the  semiempirical  codes  in  Table  1-7. 
This  will  be  followed  by  the  new  technology  developed  for  the  latest  version  of  the 
Aeroprediction  Code  (AP93).  Finally,  a  comparison  with  experiment  of  the  AP93  and  AP81 
will  be  given  for  several  missile  configurations. 
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TABLE  1-7.  CAPABILITIES  OF  EMPIRICAL  AND  SEMIEMPIRICAL  CODES 

FOR  MISSILES  (from  ref.  10) 


CONVENTIONAL 


unconventional 


ICF:  ONE  CRUCIFORM-FINNED  SECTION 
2CF:  TWO  *  '  ' 

JCF  THREE  *  *  * 

DD:  DOUBLE  DIAMETER 


CLASSICAL  BOOSTED  LUTING  AIRBREATHING 


A  A  A  A 


CODE  NAME 


INTAKES 


JCF  I  DD  ♦  JCF  I  ELLIPTIC  I  OPEN  I  CLOSED 


Coafigvraboa  type  baa  beta  coapmad 


C  Coafifarauo*  type  could  be  coaptnad 
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TABLE  1-7.  CAPABILITIES  OF  EMPIRICAL  AND  S EMIEMPIRI C AL  CODES 

FOR  MISSILES10  (CONTINUED) 


Computed 
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TABLE  1-8.  CAPABILITIES  OF  FULL  POTENTIAL  CODES  FOR  MISSILES10 


1CF:  ONE  CRUCIFORM-FINNED  SECTION 

2CF:  TWO  *  *  * 

3CF:  THREE  *  *  * 

DD:  DOUBLE  DIAMETER 

l 

CONVENTIONAL 

(drcukrcnieifira 

UNCONVENTIONAL 

CLASSICAL 

BOOSTED 

LUTING 

AOtBREATHINO 

X 

INTAKES 

CODE  NAME 

ORIGIN 

DATE 

MACH 

1CF 

2CF 

DD  +  3CF 

OPEN 

CLOSED 

SANDRAO 

SANDIA 

(USA) 

H 

198$ 

<  1 

Drag  only 

/ 

/ 

/ 

/ 

>  1 

/ 

/ 

/ 

/ 

/ 

1 

SIMP 

ROCKWELL 

(USA) 

29 

1986 

<  1 

m 

n 

/ 

/ 

/ 

/ 

>  1 

B 

B 

C 

C 

B 

C 

'/  Configuration  type  bee  been  computed  /For  no  C  Configuration  type  could  be  computed 
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TABLE  MO.  CAPABILITIES  OF  EULER  CODES  FOR  MISSILES1 


CONVENTIONAL 

(dmkfcnolm 


UNCONVENTIONAL 


1CF:  ONE  CRUCIFORM-FINNED  SECTION 
2CF:  TWO  *  *  * 

3CF:  THREE  *  *  * 

DD:  DOUBLE  DIAMETER 


CLASSICAL  I  BOOSTED  LIFTING  AIRBREATHING 


INTAKES 


CODE  NAME  |  ORIGIN  REF  N*  DATE  MACH  |  1CF  2CF  DD  +  3CF  |  ELLIPTIC  |  OPEN  |  CLOSED 


USAF 

EAGLE  ARMAMENT  43 
LAB.  (USA) 


DORNIER 


EULBMG 


EULER3D 


DORNIER 


ONERA 

•eroepelial 

(FR) 


USAF 

KODIAK  ARMAMENT  48 
Lab./Miai.  Sutc 


SANDIA 

SANDIAC  General  Electric  SO 
(USA) 


Configuration  type  hit  been  computed 


C  Configuration  type  could  be  computed 


EQUATIONS  METHOD  NUMERICAL  SCHEME  MESH  MACH 
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TABLE  1-10.  CAPABILITIES  OF  EULER  CODES  FOR  MISSILES10 

(CONTINUED) 


NOTES 

ABOUT 

SEPARATION 

(KUTTA 

Q 

1 

SSS 

SSS 

SSS 

SSS 

SSS 

SSS 

MACH 

>  1 

(2) 

* 

« 

• 

« 

• 

* 

• 

• 

• 

* 

« 

• 

• 

<  1 

(2) 

* 

* 

« 

« 

- 

• 

*■«» 

- 

- 

- 

• 

X 

3 

2 

c 

u 

e 

l 

i 

OR 

NON-STRUCT. 

(I) 

C/1 

NS 

CO 

CO 

NS 

CO 

CO 

CO 

CO 

in 

in 

CO 

CO 

NUMERICAL  SCHEME 

c 

c 

8 

\A 

1° 

SPACE 

cm 

Y 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

TIME 

- 

cm 

CM 

CM 

- 

CM 

- 

- 

- 

CM 

- 

£ 

! 

u 

C 

OR 

NON-CENT. 

0) 

NC 

NC 

u 

u 

u 

NC 

NC 

u 

u 

NC 

V 

u 

NC 

t 

2 

!•! 

2 

is 

s 

- 

E  +  1 

ui 

U1 

- 

Ul 

U 

Ul 

Ul 

Ul 

Ul 

Ul 

Ul 

METHOD 

• 

u 

2 

m 

l 

E 

i 

9 1 

j 

5  a 
E u 

» 

■ 

Is 

1 

FV 

FV 

FV 

FV 

FV 

FV 

FV 

FD 

FD 

FD 

FD 

FV 

FV 

EQUATIONS 

CONSERVATIVE 
OR  NON  CONS, 
cnau 

>  — 

U 

u 

u 

W 

o 

U 

O 

V 

U 

z 

+ 

u 

O 

u 

u 

STEADY 

OR 

UNSTEADY 

(1) 

X 

3 

3 

CO 

3 

CO 

CO 

CO 

CO 

CO 

3 

CO 

CODE 

NAME 

EAGLE 

EUFLEX 

EULBMG 

EULER3D 

EULSSM 

V 

r+s 

3 

u. 

KODIAK 

MISSILE 

MUSE 

SANDIAC 

SWINT 

WINGA2 

ZEUS 

(1)  Only  the  first  le»er(*) 

(2)  *  for  YES 

/  for  NO 

SSS  Separation  on  Smooth  Surface  with  Kutta  Condition 
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TABLE  1-11.  CAPABILITIES  OF  FULL  NAVIER-STOKES  CODES  FOR  MISSILES10 


1CF:  ONE  CRUCIFORM-FINNED  SECTION 

JCP:  TWO  *  *  • 

3CF:  THREE  •  *  ‘ 

DD:  DOUBLE  DIAMETER 

conventional 

(orator  rtow 

UNCONVENTIONAL 

CLASSICAL 

BOOSTED 

LIFTING 

AIRBREATHING 

A 

A 

X 

INTAKES 

CODE  NAME 

ORIGIN 

REF  N* 

DATE 

LAMINAR 

TURBULENT 

MACH 

1CF 

2CF 

DD  +  SCF 

ELLIPTIC 

1^21 

CLOSED 

ARC3D 

NASA 

ARC 

(USA) 

■ 

LfcT 

<  1 

■9 

C 

C 

C 

m 

IIP 

1 

& 

>  1 

| 

B 

C 

c 

m 

PSD 

NASA 

ARC 

(USA) 

1983 

L*T 

<  1 

C 

c 

c 

C 

m 

C 

m 

c 

c 

C 

B 

c 

DORNIER 

(FRG) 

1987 

<  i 

wm 

B 

c 

c 

B 

c 

>  i 

c 

c 

c 

c 

B 

c 

UWTN 

NASA 

ARC 

(USA) 

1987 

LftT 

<  i 

B 

B 

c 

c 

B 

c 

>  i 

B 

c 

c 

B 

c 

CFUDE 

NASA 

S3 

1987 

L*T 

ph 

11111 

B 

c 

c 

B 

c 

GASP 

VPI 

J4. 

J5 

1990. 

1992 

UlT 

Q 

■ 

D 

c 

c 

a 

a 

Confifurmuon  type  hee  beco  mmpuleH  C  CgeTifunbae  type  oould  be  nnmpiPri 
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TABLE  1-12.  CAPABILITIES  OF  PARABOLIZED  NAVIER-STOKES  CODES 

FOR  MISSILES10 


ICF  ONE  CRUCIFORM- FINNED  SECTION 

2CF:  TWO  *  ■  * 

JCF  THREE  •  '  ' 

DEE  DOUBLE  DIAMETER 

CONVENTIONAL 

UNCONVENTIONAL 

CLASSICAL 

BOOSTED 

LIFTING 

AIRBREATHING 

A 

X 

INTAKES 

CODE  NAME 

ORIGIN 

REF  N* 

DATE 

LAMINAR 

TURBULENT 

MACH 

ICF 

2CF 

DD  +  JCF 

ELLIPTIC 

OPEN 

CLOSED 

PNS 

(basic  venian) 

$6 

57 

1979 

1984 

LAT 

>  1 

■ 

a 

C 

*  -  ~ 

D 

n 

PNSFVM 

DORN1ER 

(FRO) 

1916 

L 

>  1 

H 

B 

C 

a 

n 

CoaTifurmUoo  type  bam  been  computed  C  CooftfurtlMMi  type  could  be  computed 
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2.0  CONVENTIONAL  APPROXIMATE  AERODYNAMIC  METHODS 


This  section  of  the  report  will  review  some  of  the  more  important  approximate 
aerodynamic  methods  that  have  proved  quite  useful  in  the  development  of  semiempirical 
codes.  Time  and  space  will  not  permit  derivation  of  the  methods  from  first  principles. 
However,  appropriate  references  will  be  given  for  the  interested  reader.  The  approach  taken 
here,  in  the  presentation  of  the  material,  will  be  to  mention  the  assumptions  inherent  in  each 
method,  relevant  equations,  and  possibly  show  an  example  or  two  as  may  be  warranted. 


2. 1  HYBRID  THEORY  OF  VAN  DYKE  (HTVD)58 


The  Hybrid  Theory  of  Van  Dyke5*  combines  a  second-order  solution  to  the  potential 
equation  with  a  first-order  crossflow  solution  first  espoused  by  Tsien59.  The  advantage  of 
this  method  is  that  it  gives  second-order  accuracy  in  the  axial  direction  where  first-order 
accuracy  is  generally  unacceptable  for  drag  computations.  On  the  other  hand,  first-order 
accuracy  in  the  crossflow  plane  is  typically  acceptable  for  normal  force  and  center  of 
pressure  computations.  The  fundamental  reason  for  this  is  that  perturbations  in  the  flow,  due 
to  the  presence  of  a  body,  have  more  impact  in  the  axial  as  opposed  to  the  normal  force 
direction.  Hence,  to  get  axial  force  accuracy  compatible  with  a  goal  of  ±  10  percent  require 
second-order  methods,  whereas  ±  10  percent  accuracy  on  CN  can  be  obtained  with  first-order 
methods  in  many  cases. 

As  already  mentioned,  the  Hybrid  theory  comes  from  the  potential  equation  of  fluid 
mechanics.  It  is  limited  to  supersonic  flow  (we  have  used  this  method  down  to  M„  =  1.2) 
where  the  assumption  of  ise  itropic  flow  (shock  waves  are  weak)  can  be  made.  This  typically 
limits  the  upper  Mach  number  range  to  about  M.  =  2.0  to  3.0,  depending  on  the  body 
shape.  Also,  the  slope  of  the  body  surface  must  be  less  than  the  Mach  Angle.  The  Tsien 
solution,  or  crossflow  part  of  the  solution,  comes  from  the  linearized  perturbation  equation. 
On  the  other  hand,  the  second-order  solution  to  the  axial  flow  is  found  by  obtaining  a 
particular  solution  to  a  reduced  version  of  the  full  potential  equation.  This  is  the  key  to  the 
accuracy  improvement  afforded  by  Van  Dykes  solution  in  that  some  of  the  nonlinearity 
inherent  in  the  axial  flow  problem  is  brought  into  the  solution  by  this  process.  The  beauty  of 
the  Van  Dyke  method  is  that  this  particular  second-order  solution  is  given  entirely  in  terms 
of  the  first-order  solution.  That  is,  one  simply  solves  the  first-order  perturbation  solution  for 
the  axial  flow  and  then  solves  an  algebraic  equation  for  the  second-order  solution  where  the 
boundary  condition  at  the  body  is  satisfied. 


2-1 


NSWCDD/TR-93/55 1 


In  equation  form,  the  general  first-order  perturbation  problem  is  5*: 

$rr  +  $r/r  +  ®M/x2  -  (A£?  -  DGjo,  =  0  (1) 

with  boundary  conditions  that  do  not  allow  any  upstream  disturbances: 

$  (0,r,$)  =  Qx  (0,r,4>)  =  0  (la) 

and  that  require  the  flow  to  be  tangent  to  the  body  surface: 

0r  (x,rb,e)  +  sina  cose  *  ■—  [cosa  +  (x,  zb,z)]  (lb) 


The  subscripts  in  Equation  (1)  indicate  partial  derivatives.  The  solution  to  Equation  (1)  is 
satisfied  identically  by: 

<t>  {x,r,o)  -  Y1(x,r)  cosa  +  Ci  ( x,r )  sina  cose  (2) 

The  first  term  of  Equation  (2)  is  the  first-order  axial  solution,  and  the  second  term  is  the 
first-order  crossflow  solution.  Since  the  equation  is  linear,  these  two  solutions  can  be  found 
independently,  and  then  added  together.  The  axial  solution,  Y,  (x,  r),  for  a  general  body  is 
found  by  placing  a  series  of  sources  and  sinks  along  the  x  axis  and  satisfying  the  boundary 
conditions  at  each  point.  The  crossflow  solution,  f,(x,  y),  is  found  by  placing  a  series  of 
doublets  along  the  axis,  again  satisfying  the  boundary  conditions. 

The  particular  second-order  solution  that  Van  Dyke  found  for  the  reduced  full 
potential  equation  is 


Y2  =  id  [Y^  (Yx  +  tfrYlr)  -  (■£)  Y3lr]  where  N  =  (X_±_i)  fjjL  (3) 

Second-order  axial  velocity  components  and  Y2r  are  also  defined  in  terms  solely  of  the 
first-order  solution  Y,(x,r). 
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Once  the  second-order  axial  perturbation  velocity  components  ^2r  are  computed, 
along  with  the  first-order  crossflow  components  fu  and  flr,  the  total  perturbation  velocities 
are  then: 

=  (cosa)  (l+W2x)  +  (sina  cose)  Cut  (4a) 


-~r  =  cosa  (Y2r)  +  (sina  cose)  (l  +  £lr) 


—  *  -  (sina  sine)  (l  +  -ii) 
Vm  r 


:i  ♦  ii, 


The  pressure  coefficient  at  each  body  station  is  then: 


(*,»)  -4:  1*  jg  1  -  V‘-+-W  (5) 


Finally  the  force  coefficients  are: 


-  2  r‘  r 

nr?  Jo  Jo 


Cp  ( x ,  o)  ~L  do  dx 


cN  =--2—  f  f  Cp  (x,  e)  cos  (o)r  de  dx 
nr t  Jo  Jo 


cm  =  — “~7  f  f  Cp  (x,  e)  cos  (e)x  r  de  dx 
nr:  Jo  Jo 


and  the  center  of  pressure  in  calibers  from  the  nose  is 


*C,  " 


It  should  be  pointed  out  that  in  the  actual  numerical  integration  of  Equations  (6),  (7),  and 
(8)  the  integration  must  be  carried  out  in  segments  of  the  body  between  each  discontinuity 
due  to  the  discontinuous  pressure  distribution. 
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Also,  the  hybrid  theory  of  Van  Dyke  is  limited  to  pointed  bodies  of  revolution. 
Bluntness  will  be  considered  later. 


2.2  SECOND-ORDER-SHOCK-EXPANSION  THEORY  (SOSET) 


First-order  Expansion  Theory  was  first  proposed  by  Eggers  et  al.  for  bodies  of 
revolution  flying  at  high  supersonic  speeds.61"63  Basically,  the  Shock-expansion  Theory 
computes  the  flow  parameters  at  the  leading  edge  of  a  two-dimensional  (2-D)  surface  with 
the  oblique  shock  wave  relations  and  with  the  solution  for  a  cone  at  the  tip  of  a  three- 
dimensional  (3-D)  body.  Standard  Prandtl-Meyer  Expansion  (PME)  is  then  applied  along  the 
surface  behind  the  leading  edge  or  tip  solution  to  get  the  complete  pressure  distribution  over 
the  body  surface.  Referring  to  Figure  2-1,  this  theory  inherently  assumes  that  the  expansion 
waves  created  by  the  change  in  curvature  around  the  body  are  entirely  absorbed  by  the  shock 
and  do  not  reflect  back  to  the  body  surface.  Since  the  theory  assumes  constant  pressure 
along  one  of  the  conical  tangent  elements  of  the  surface,  fairly  slender  surfaces  must  be 
assumed  or  many  points  along  the  surface  assumed  to  obtain  a  fairly  accurate  pressure 
distribution.  Another  way  of  stating  this  is  to  minimize  the  strength  of  the  disturbance 
created  by  Mach  waves  emanating  from  the  expansion  comer  and  intersecting  the  shock,  the 
degree  of  turn  should  be  small. 


FIGURE  2-1.  APPROXIMATION  OF  TRUE  BODY  BY  ONE  COMPOSED 
OF  STRAIGHT  LINE  SEGMENTS  TANGENT  TO  SURFACE 
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Syvertson  (et  al.)  extended  the  generalized  Shock-expansion  Theory  on  pointed  bodies 
and  sharp  airfoils  to  what  he  called  a  second-order  theory.80  He  defined  the  pressure  along  a 
conical  frustum  by 

P  *  Pc  ~  (Pc  “  P2)  (10) 

instead  of  a  constant  on  each  segment  as  was  the  case  in  the  generalized  theory.  Here  Pc  is 
the  pressure  on  a  cone  with  the  given  cone  half  angle  equal  to  the  slope  of  the  conical 
segment  with  respect  to  the  axis  of  symmetry.  pj  is  the  pressure  just  aft  of  a  conical  segment 
(see  Figure  1)  which  is  calculated  from  a  Prandt  Meyer  Expansion  (PME)  of  the  flow  around 
a  comer  (as  shown  in  Figure  2-2,  going  from  points  1  and  3  to  points  2  or  4,  for  example). 


Also 


FIGURE  2-2.  FLOW  ABOUT  A  FRUSTUM  ELEMENT 


►  x 


(10a) 


Thus,  examining  p  from  Equation  (10),  it  can  be  seen,  for  example,  on  the  ffustrum  element 
in  Figure  2-2  that  the  pressure  varies  from  the  pressure  ot  the  generalized  theory  at  point  2 
to  that  of  a  cone  of  angle  62  and  Mach  number  M2  as  s  gets  large.  Syvertson  and  Dennis 
approximated  the  pressure  gradient  as80 


(M 
o2  las), 


(id 
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where 


=  YPl,2^1.2 

2  (Afi.a  -  1) 


Q 


1 , 2  = 


1  + 


111 

2 


y*1 

2(Y-1> 


Finally,  for  negative  angles  such  as  would  occur  on  a  boattailed  configuration,  pc  was 
replaced  by  p*.  No  discussion  was  given  for  blunt  bodies.  It  should  be  noted  that  if  tj  of 
Equation  (10)  becomes  negative,  the  SOSET  reverts  to  the  generalized  or  first-order  Shock- 
expansion  Theory.  This  is  because  Equation  (10)  will  not  give  the  correct  asymptotic  cone 
solution  for  negative  values  of  t\. 

Experience  has  shown  that  SOSET  gives  very  good  pressure  distributions  for  low  to 
moderate  angles  of  attack  and  at  Ma  2:  2.  As  Mach  numbers  decrease  below  about  2.5,  the 
SOSET  becomes  increasingly  inaccurate  until  about  M*  =  1.5,  where  the  accuracy  is 
generally  unacceptable.  This  applicable  Mach  number  range  is  very  complimentary  to  the 
Hybrid  Theory  of  Van  Dyke  where  the  accuracy  is  best  between  1.2  5  M.  5  2.5. 


2.3  MODIFIED  NEWTONIAN  THEORY  (MNT) 


Newtonian  Impact  Theory  assumes  that,  in  the  limit  of  high  Mach  number,  the  shock 
lies  on  the  body.  This  means  that  the  disturbed  flow  field  lies  in  an  infinitely-thin  layer 
between  the  shock  and  body.  Applying  the  laws  of  conservation  of  mass  and  momentum 
across  the  shock  yields  the  result  that  density  behind  the  shock  approaches  infinite  values  and 
the  ratio  of  specific  heats  approaches  unity.  The  pressure  coefficient  on  the  surface 
becomes64 


Cp  =  2sin26  e(J  (12) 

where  8^  is  the  angle  between  the  velocity  vector  and  a  tangent  to  the  body  at  the  point  in 
question  (see  Figure  2-3).  8^  is  defined  by: 

sin  (6S7)  =  sine  sina  -  sina  cose  cos0  (13) 
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Lees12  noted  that  a  much  more  accurate  prediction  of  pressure  on  the  blunt-nose  body 
could  be  obtained  by  replacing  the  constant  "2“  in  Equation  (12)  with  the  stagnation  pressure 
coefficient  CPe.  CPo  can  be  calculated  from: 


'Po 


2 

Y  Mi 


’  (Y  +  1)a£ 

7^~i 

Y  +  1  1 

i 

T  -  1  -  1 

2 

2yt£  -  (Y  -  1). 

X 

(14) 


MNT  is  thus  defined  by: 

CP  =  CJ>„  Sin26w  (15) 

Equation  (IS)  allows  the  calculation  of  the  pressure  coefficient  all  along  the  blunt  surface  of 
a  missile  nose  or  wing  leading  edge  for  a  perfect  gas  where  CPo  is  given  by  Equation  (14) 
and  sin  8^  from  Equation  (13). 


z 


z 


WHERE:  r2  -  y2  ♦  z* 


FIGURE  2-3.  NOMENCLATURE  USED  FOR  DETERMINATION  OF  ANGLE  8 „ 

Experience  has  shown  that  the  MNT  gives  very  acceptable  estimates  of  pressure 
coefficient  on  the  blunt  portion  of  a  nose  or  leading  edge,  even  at  Mach  numbers  where  the 
assumptions  of  Newtonian  Impact  Theory  are  violated. 
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2.4  HYBRID  THEORY  OF  VAN  DYKE  COMBINED  WITH  MODIFIED 
NEWTONIAN  THEORY  (HTVD/MNT)1 


As  noted  in  the  discussion  on  the  Hybrid  Theory,  it  is  limited  to  conditions  where  the 
body  slope  is  less  than  the  local  Mach  angle.  This  means  it  is  not  applicable  in  the  nose 
region  of  a  blunt  missile.  On  the  other  hand,  MNT  gives  very  acceptable  estimates  of 
pressure  coefficients  in  the  nose  region,  even  for  low  supersonic  Mach  numbers  where  the 
assumptions,  inherent  in  the  Newtonian  Impact  Theory,  are  violated.  Moore1  was  the  first  to 
recognize  the  possibility  of  combining  these  two  theories.  The  key  to  the  successful 
combination  was  in  the  starting  solution.  At  low  supersonic  Mach  numbers,  the  pressure 
overexpands  on  a  blunt  nose  tip  as  it  proceeds  around  the  blunt  portion  from  the  stagnation 
point  to  the  given  portion  of  the  nose.  In  order  to  capture  this  overexpansion,  Moore  found 
that  it  was  necessary  to  start  the  HTVD  near  its  maximum  acceptable  slope  and  allow  the 
pressure  to  expand  around  the  surface.1  Simultaneously,  the  MNT  was  started  at  the 
stagnation  point  and  allowed  to  expand  until  the  pressure  coefficients  of  the  MNT  and  the 
HTVD  were  equal.  This  was  defined  as  the  Match  point.  Upstream  of  the  Match  point, 
MNT  was  used  in  the  force  and  moment  calculations,  whereas  downstream,  HTVD  was 
used.  Figure  2-4  is  an  illustration  of  the  boundaries  of  perturbation  and  Newtonian  theories. 
Figures  2-5  and  2-6  illustrate  the  capability  of  this  theory  to  accurately  predict  pressure 
coefficients  on  a  35  percent  blunt  cone  of  11.5°  half  angle  at  a  =  8°  and  at  Mx  =  1.5  and 
2.96.  Note  the  excellent  agreement  of  the  combined  theory  all  along  the  surface  at  M,  = 
1.5.  Particularly  impressive  is  its  ability  to  capture  the  overexpansion  region  around  x  = 

0.1  to  x  =  0.4.  Also,  note  that  SOSET  gives  fairly  poor  estimates  at  M.  =  1.5.  On  the 
other  hand,  at  M*  =  2.96,  the  HTVD/MNT  is  no  better  (and  mavbe  slightly  worsen  than  the 
SOSET/MNT,  which  will  be  discussed  next. 

To  the  authors  knowledge,  the  HTVD/MNT  remains  the  only  accurate  engineering 
method  to  estimate  low  supersonic  Mach  number  aerodynamics  for  blunt  and  sharp  tip  bodies 
of  revolution.  Attempts  were  made  to  extend  the  SOSET/MNT  down  to  the  low  supersonic 
Mach  number  range,  but  without  success. 
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FIGURE  2-4.  BOUNDARIES  OF  PERTURBATION  AND  NEWTONIAN  THEORY 
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FIGURE  2-5.  COMPARISON  OF  THEORY  AND  EXPERIMENT  FOR 

BLUNTED  CONE 


—  =  0.35,  K.  =  1.5,  a  =  8°, ec=li.5° 
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FIGURE  2-6.  COMPARISON  OF  THEORY  AND  EXPERIMENT  FOR 

BLUNTED  CONE 
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2.5  SECOND-ORDER-SHOCK-EXPANSION  THEORY  COMBINED  WITH 

MODIFIED  NEWTONIAN  THEORY  (SOSET/MNT)65-68 

Jackson  et  al.65  combined  SOSET  with  MNT  to  treat  blunt-nosed  configurations  with 
or  without  flares.  Jackson  et  al.65,  like  Syvertson  and  Dennis60,  assumed  that  the  lifting 
properties  could  be  predicted  by  assuming  that  the  original  body  is  made  up  of  several 
equivalent  bodies  of  revolution  represented  by  the  various  meridians  (see  Figure  2-7).  They 
assumed  the  match  point  between  the  MNT  and  second-order  shock  pressure  prediction  to  be 
the  angle  that  corresponds  to  shock  detachment  on  a  wedge  with  the  given  freestream  Mach 
number. 

De  Jamette  et  al.66*8  made  significant  improvements  to  the  work  of  Jackson  et  al.65 
and  Syvertson.60  These  new  improvements  included  the  following: 

1.  An  exact  (as  opposed  to  an  approximate)  expression  for  the  pressure  gradient 
downstream  of  a  comer. 

2.  A  new  expression  for  pointed-cone  pressures  at  angle  of  attack  which  improves 
the  initial  pressure  prediction  over  that  of  tangent  cone  theory. 

3.  A  new  technique  for  calculating  pressures  on  bodies  at  incidence. 

The  pressure  computations  at  angle  of  attack,  showed  improvement  over  the  method 
of  Jackson65,  De  Jamette,  el  al.23  derived  a  new  expression  for  pointed-cone  pressure  at 
a  >  0  by  combining  Slender  Body  Theory,  Newtonian  Theory,  and  an  approximate 
expression  for  Cp^  to  give: 


Cp(a,0,z,M)  =  CPmwo  +  ACp 


(16a) 


where 


ACp=  -sin2asin20cos$+sin2acos20  (l-tan20)  -|2  +  -^jsin2® 


(16b) 


Cp  =  sin20 
*■•0  c 


1+ 


(v  +  l)*2  +  2 
(y  -  1)  K2  +  2 


to 


/T  1 1  ♦  JL 

\  2  K2 


(16c) 


and 

K2  =  (id  -  1)  sin20c 
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FIGURE  2-7.  TYPICAL  EQUIVALENT  BODY  SHAPES  USED  FOR  COMPUTING 
LIFTING  PROPERTIES  WITH  SECOND-ORDER  SHOCK  EXPANSION  THEORY 
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Note  also,  that  while  Equation  16  was  strictly  defined  for  pointed  cone  pressures  at  angle  of 
attack,  it  could  also  be  used  in  a  Tangent  cone  sense  to  obtain  pressures  at  any  point  on  a 
body  surface.  De  Jamette  actually  used  loading  functions  to  obtain  body  alone  lift 
properties,  however6®. 

Figure  2-8  presents  results  of  De  Jamette  et  al  “  compared  to  experiment.  The  case 
chosen  is  the  same  case  shown  in  Figure  2-6,  except  here,  the  method  of  De  Jamette  et  al 68 
is  used  versus  Jackson  et  al63  in  Figure  2-6.  While  the  two  theories  are  close,  in  comparing 
Figure  2-8  with  Figure  2-6,  it  is  seen  that  the  theory  of  De  Jamette  et  al65  does  show 
improvement  in  pressure  prediction  and  therefore  forces  and  moments  as  well. 


2.6  ALLEN-PERKINS  VISCOUS  CROSSFLOW  THEORY69 


A  fairly  simple,  yet  quite  powerful,  method  for  computing  body-alone  nonlinear 
aerodynamics  was  introduced  by  Allen -Perkins69.  Allen  reasoned  that  the  total  force  on  an 
inclined  body  of  revolution  is  equal  to  the  potential  term  discussed  previously  plus  a  cross 
flow  term.  This  term  is  based  on  the  drag  force  experienced  by  an  element  of  a  circular 
cylinder  of  the  same  diameter  in  a  stream  moving  at  the  cross  component  of  the  stream 
velocity,  V.  sin  a.  This  crossflow  term  is  primarily  created  by  the  viscous  effects  of  the 
fluid  as  it  flows  around  the  body,  often  separating  and  creating  a  nonlinear  normal  force 
coefficient.  In  equation  form,  the  so  called  viscous  crossflow  theory  is: 


(17) 


Here  ij  is  the  drag  proportionality  factor  or  crossflow  drag  of  a  cylinder  of  finite  length  to 
one  of  infinite  length.  Cd<_  is  the  crossflow  drag  coefficient.  Also,  the  crossflow  theory 
assumes  the  center  of  pressure  of  the  nonlinear  term  is  at  the  centroid  of  the  planform  area. 
Generally,  the  total  center  of  pressure  is  a  weighted  average  of  the  linear  and  non  linear 
components  of  normal  force.  That  is 


V  _  (Xc CHaL  +  (Xcp>  L  C«L 
Acp  F" 


(18) 


The  pitching  moment  about  a  given  point  Xq  is  then 

=  -CN(Xcp  "  V  (19) 

The  original  work  of  Allen  did  not  include  compressibility  effects  in  n  but  Reynolds  number 
effects  were  shown  in  Cde  at  low  crossflow  Mach  numbers. 
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FIGURE  2-8.  PRESSURE  DISTRIBUTION  ON  A  BLUNTED  CONE 
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2.7  VAN  DRIEST  U  METHOD  FOR  SKIN  FRICTION  DRAG70 


Another  powerful,  yet  simple,  method  for  performing  computation,  is  the  Van  Driest 
II  method  for  computing  skin-friction  drag.  This  method,  as  derived,  is  based  on  two 
dimensional  turbulent  boundary  layer  flow.  Strictly  speaking,  it  is  only  applicable  to  regions 
of  flow  on  the  lifting  surfaces  where  the  flow  is  turbulent,  two  dimensional,  and  the  viscous 
region  is  primarily  confined  to  a  thin  layer  near  the  surface  (boundary  layer).  In  practice, 
however,  it  has  been  applied  to  two  and  three  dimensional  surfaces  with  success. 

The  turbulent  mean  skin-friction  coefficient  according  to  Van  Driest70  is: 

— 7— 172  =  ( VrJ  1/2  (sin^q+sin-1^)  *log10  ( R0mCfJ 

(20) 

-  (^)lo9io(Vrj 


where 


2  A*-B 


(B2+4A2)1/2 


B 

(B2+4A2)1/2 


and 


A  = 


|l/2 


l+(y-l)/2Ad 


The  variable  n  of  Equation  (20)  is  the  power  in  the  power  viscosity  law: 


JL  = 

H- 


(21) 
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The  freestream  Reynolds  number  and  adiabatic  wall  temperature  are  given  by: 


P,  V.  t 
p- 


(22) 


~  =  1  ♦  0.9  X_=-A  (23) 

Tm  2 

Equations  (20)  through  (22)  allow  the  calculation  of  the  mean  turbulent  skin-friction  over  the 
entire  body  or  wing  area.  The  skin-friction  axial  force  coefficient  on  each  component  is 
then: 


Cl,  -  cf_  ^ 

Az»f 


(23) 


where  AweJ  is  the  surface  area  of  the  component  in  question. 

For  most  flows,  a  portion  of  the  flow  is  laminar.  An  approximation  to  the  mean 
skin-  friction  coefficient  for  laminar  flow  can  be  obtained  from70 


1.328 


(24) 


Here  the  Reynolds  number  is  based  on  the  distance  where  transition  occurs  rather  than  the 
reference  length,  as  was  the  case  for  Equation  (22). 

The  point  where  transition  occurs  is  dependent  on  many  factors.  Experience  has 
shown,  for  flight  vehicles,  a  transition  Reynolds  number  of  1  x  106  for  the  body  and  0.5  x 
10®  for  the  wings  gives  acceptable  numbers.  For  wind  tunnel  models  without  a  trip,  a 
transition  Reynolds  number  of  3  to  5  million  is  more  reasonable  due  to  a  smooth  surface.  If 
a  boundary  layer  trip  is  used,  the  entire  configuration  component  should  have  turbulent  flow. 
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2.8  LIFTING  SURFACE  THEORY71 

Lifting  Surface  Theory  refers  to  the  solution  of  the  flow  over  a  three  dimensional 
wing  where  the  distribution  of  pressure  is  allowed  to  vary  in  both  the  spanwise  and 
chordwise  direction.  The  fundamental  equation  is  the  three  dimensional  pemirbation 
equation,  here  written  in  rectangular  coordinates,  as: 

(1  -Mi)  -  0  (25) 


The  Flow  tangency  boundary  condition  requires: 


_  dzu 
■ 

dzx 

= 


at  z  *  0* 
at  z  -  0 


>  for  {x,  y)  on  S 


(25a) 


If  the  wing  thickness  is  neglected  and  we  limit  ourselves  to  missiles,  then  wing  chamber  can 
also  be  neglected.  Then  the  boundary  conditions  in  Equation  (25a)  become: 

*  -o  (25b) 


for  both  the  upper  and  lower  surfaces. 

In  addition  to  this  boundary  condition,  the  Kutta  condition  (which  requires  the  velocity  on  the 
upper  and  lower  surfaces  at  the  trailing  edge  to  be  equal)  is  also  imposed  for  subsonic  flow. 

The  assumptions  involved  in  the  Lifting  Surface  Theory,  as  applied  to  most  missile 
configurations,  are  therefore  small  perturbations  in  the  flow  due  to  the  presence  of  the  wing 
and  the  thickness  and  chamber  effects  are  zero  or  small  compared  to  angle  of  attack  effects. 
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Equation  (25)  may  be  simplified  somewhat  by  using  Prandtl-Glauert  rule  (72)  to 
relate  the  compressible  subsonic  normal  force  or  pitching  moment  to  the  incompressible  case. 
That  is: 


{CN) 


K..A *.« 


(  ^  0 ,  AK, « 

sjl  -  Ml 


( CM) 


H..AS.U 


( Cg}  o,As,m 
-  M? 


(26) 


Using  the  above  relations,  the  normal  force  and  pitching  moment  on  a  given  wing  at  any 
subsonic  Mach  number  may  be  found  by  calculating  the  aerodynamics  of  the  same  wing  at 
zero  Mach  number.  Figure  2-9  is  a  representation  of  the  wing  planform  parameters. 


For  M„  =  0,  Equation  (25)  reduces  to  La  Places  equation 


V2®  =  0 


(27) 


with  boundary  condition  (25b). 


FIGURE  2-9.  WING  PLANFORM  PARAMETERS 
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There  are  many  methods  to  solve  Equation  (27).  The  one  used  here  is  that  of 
Chadwick  et  al.71  which  closely  follows  Ashley  et  al.72.  The  velocity  potential  $  is  given 
by: 


Q{x,y,  z) 


-L  1 

■fACp(Xl,yQ  z 

x-x, 

1+  1 

8*  { 

r  ( y-y i)2+*2 

^(jr-x1)2+(y-y1)2+z2 

dx1dy1(28) 


Here,  x„  y,  are  coordinates  of  an  element  of  the  lifting  surface  that  has  a  differential 
pressure  coefficient  of  ACP  between  the  lower  and  upper  surfaces  at  this  point  (x,,  y,).  It  is 
required  to  determine  the  pressure  loading  over  the  entire  surface.  Following  Chadwick71, 
Equation  28  is  first  differentiated  with  respect  to  z  and  the  limit  as  z  -*  0  taken.  The  result 
is  then  equated  to  the  boundary  condition,  Equation  (25b)  to  obtain: 


f  ACJp(x1,y1) 

i+  x  Xl 

/  (y-yj 2 

Jix-xJ**  (y-yj)2 

dx1dy1  (29) 


The  cross  on  the  y,  integral  indicates  a  singularity  at  y  =  y,,  in  which  case  Manglers 
principal-value  technique72  can  be  applied.  The  details  of  the  solution  of  the  integral 
Equation  (29)  for  ACP  (x,y)  will  not  be  repeated  here  as  they  are  given  in  detail  in  many 
references  (see  for  example,  Chadwick.71)  Worthy  of  note,  however,  is  the  fact  that 
Equation  (29)  is  an  integral  equation  for  which  the  wing  loading  AC,,  is  to  be  found  as  a 
linear  function  of  angle  of  attack.  This  wing  Loading  is  first  approximated  by  a  series 
expansion  with  a  set  of  unknown  coefficients  of  number  equal  to  the  number  of  surface 
elements  on  the  wing  planform.  That  allows  each  ACP  to  be  influenced  by  all  other 
elements  of  the  wing.  The  unknown  coefficients  in  each  ACP  series  are  found  by  solution  of 
an  inverse  matrix.  ACP  (x,y)  is  then  calculated. 


Once  the  span  loading  ACP  (x,y)  is  known  over  the  entire  wing  surface,  the  normal 
force  at  a  given  spanwise  location  is: 


TE 

XLB 


The  total  normal  force  for  the  entire  wing  is: 

c  =  r 

^ref  j  I 


A  Cpdx 


>b/2 

ref  J  0 


cC„dy 


(30) 


(31) 


The  pitching  moment  of  a  given  airfoil  section,  about  the  point  where  the  wing  leading  edge 
intersects  the  body,  is  then  (positive  leading  edge  up): 
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= 


1  f XTE 

C^z0f  J  xLE 


xA  Cpdx 


(32) 


The  total  pitching  moment  become*' 


Cu  = 


2  f*/2 


Sr»f  J  o 


/: 


cCmdy 


(33) 


If  it  is  desired  to  calculate  the  pitching  moment  about  some  other  reference  point,  then 


C  ♦  C  ^ 


(34) 


z»f 


where  Xo  is  the  distance  from  the  reference  point  to  the  juncture  of  the  wing  leading  edge 
with  the  body.  The  center  of  pressure  of  an  airfoil  section  is: 


or  of  the  entire  wing 


(35) 


P6) 


Finally,  the  span  wise  center  of  pressure  of  a  wing  semispan  is: 


(37) 


Equations  (30),  (31),  (32),  (33),  and  (37)  can  be  solved  by  numerical  quadrature,  such  as 
Simpson’s  rule,  with  special  attention  given  to  the  leading  edge  singularity. 


It  should  also  be  mentioned  that  if  one  is  interested  in  dynamic  derivatives73,  these 
aerodynamics  can  be  obtained  by  a  modification  to  the  boundary  condition,  Equation  (25a). 
That  is,  for  rolling  and  pitching  motions,  the  angle  of  attack  in  Equation  (25a)  is  replaced 
by: 


a  (x,  y) 


.♦a**  ^ 

o  y  y 


(38) 
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Equation  (27)  is  a  linear  partial  differential  equation  so  that  solutions  can  be  combined 
together  in  a  linear  fashion.  This  means,  for  roll  damping,  simply  set  aD  =  q  =  0  and  the 
boundary  condition  is 

ttU.y)  =  (38a) 


Likewise,  for  pitching  damping,  a0  =  p  =  0  and 

.  .  q{x-xraf) 

a(x,y)  =  ■-  vref/  (38b) 


2.9  THREE  DIMENSIONAL  THIN  WING  THEORY72 


Three  Dimensional  Thin  Wing  Theory  (TDTWT)  is  quite  similar  to  lifting  surface 
theory  (LST)  in  the  sense  the  same  perturbation  Equation  (25)  is  used.  The  only  difference 
is  that  TDTWT  is  normally  used  to  represent  the  supersonic  flow  solutions  of  Equation  (25) 
versus  LST  the  subsonic  solutions.  Since,  for  supersonic  flow,  solutions  to  Equation  (25) 
are  hyperbolic  versus  elliptic  for  the  subsonic  case,  they  generally  are  easier  to  obtain.  This 
is  because  no  upstream  influence  is  felt  by  a  disturbance  at  a  given  point  on  the  wing 
surface.  In  contrast,  the  subsonic  solutions  required  a  matrix  inversion  at  each  wing  element 
to  determine  the  unknown  coefficients  used  to  determine  the  pressure  differential  from  lower 
to  upper  surfaces.  On  the  other  hand,  the  assumptions  of  TDTWT  are  the  same  as  for  LST. 
They  both  assume  small  perturbations  in  an  isentropic  flow.  The  isentropic  flow  assumption 
means  no  shock  waves  are  allowed. 

In  contrast  to  the  body  solutions  generated  by  Van  Dyke,  adequate  wing  solutions  can 
be  obtained  at  higher  Mach  numbers.  This  is  because  of  the  low  slopes  present  on  most 
wing  planforms  (thickness  is  generally  very  small),  the  wing  frontal  area  is  generally  less 
than  10  percent  of  the  body  frontal  area,  and  in  the  region  of  leading  edge  bluntness,  where 
perturbation  theory  is  invalid,  modified  Newtonian  Theory  is  used  for  wave  drag  calculation. 

The  most  general  boundary  conditions  for  Equation  (25)  in  supersonic  flow  are  the 
flow  tangency  condition  specified  by 


v(x,  y)  _ 

Vm 


- 


dF 

Tx 


■m 


+  o  + 


*.y 


EL 

v_ 


q{x-x[0f) 

vm 


at  (39) 


and  the  perturbation  velocities  must  vanish  upstream  from  the  point  where  the  disturbance 
originates.  Mathematically,  this  can  be  stated  in  the  form 
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u(o~ty,z)  =  v(o~,y,z)  =  w{o~,y,z)  =  0 


(40) 


Since  Equation  (25)  is  linear,  individual  solutions  can  be  added  together.  This  allows 
individual  treatment  of  the  Equation  (39)  boundary  condition  for  drag,  lift,  roll  and  pitch 
damping  computations.  For  wave  drag  calculations,  only  the  first  term  of  Equation  (39)  is 
retained  and  the  other  terms  are  set  to  zero.  For  lift  calculations,  the  angle  of  attack  a  is 
retained  and  the  other  terms  set  to  zero.  For  roll  damping,  the  third  term  of  Equation  (29) 
is  retained  and  the  other  terms  set  to  zero.  For  pitching  rate,  the  q  term  of  Equation  (39)  is 
retained  and  the  other  terms  set  to  zero.  Finally,  for  a  constant  vertical  acceleration,  the  last 
term  is  retained  and  the  other  four  terms  set  to  zero.  Pitch  damping  moment,  CM^  +  , 

normally  refers  to  the  sum  of  the  terms  due  to  a  constant  pitch  rate  and  constant  vertical 
acceleration. 

The  solution  to  Equation  (25),  using  the  first  term  of  Equation  (39)  as  the  boundary 
condition,  will  give  the  axial  force  coefficient  of  a  sharp  wing.  If  the  leading  edge  is  blunt, 
MNT  is  used  in  conjunction  with  perturbation  theory.  The  general  solution  to  Equation  (25) 


«  (x,y,  0)  =  f  f  — 


dx1dy1 

2(y-y i)2 


(41) 


The  pressure  coefficient  at  any  point  on  the  wing  surface  is 


Cp  =  -2  Qx{x,y,o) 


(42) 


The  perturbation  velocity  at  a  given  point  p,  is  dependent  on  the  location  of  the  point 
with  respect  to  the  line  of  sources  and  sinks  which  generates  the  wing  leading  edge  or  other 
discontinuity  and  whether  this  point  is  in  a  subsonic  or  supersonic  flow  region.  For 
example,  referring  to  Figure  2-10A,  if  point  P  is  at  P,,  and  the  wing  generator  is  a  subsonic 
source  or  sink  line  (SOSL),  then 
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(43) 
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FIGURE  2-10A.  TRIANGULAR  SURFACE  SYMMETRIC  ABOUT  X  AXIS 

FOR  SUBSONIC  SOSL 


where  w  is  determined  from  the 

y  =  y„i): 


boundary  condition  and  is  (for  the  airfoil  section  at 


"(Xpi'Yei) 
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In  Equation  (43),  the  definitions 


k  =  tan  A 


(43a) 


have  been  used.  If  P  =  P2,  the  induced  velocity  at  P2  due  to  a  given  SOSL  is: 


<X> 
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2  w(xp2,yp2) 
itPv/rp^1! 


cosh'1 
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l-O2 


(44) 


At  the  wing  tip,  there  is  an  additional  disturbance  within  the  Mach  line  emanating  from  the 
tip  leading  edge  (Figure  2-10B).  The  induced  velocity  in  this  region,  P  =  P3  is: 


d> 
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w^p3,yp3) 

*PvV~l 


cosh'1 


_a!d£L 

T)  (|of+l) 


(45) 


The  absolute  value  of  a  is  taken  because  a  is  actually  negative  for  the  point  P3  .  The  induced 
velocity  at  any  point,  say  P  =  P4,  outside  of  the  Mach  lines  emanating  from  the  beginning  of 
the  SOSL  is  zero  since  this  point  is  out  of  the  zone  of  influence. 


MACH  LINE 


FIGURE  2-10B.  WING  TIP  EFFECTS  FOR  SUBSONIC  SOSL 
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If  the  wing  generator  is  supersonic,  the  Mach  lines  from  point  0  in  Figure  2-11 A  lie 
behind  the  SOSL.  If  in  Figure  2-11  A,  P  =  Pj,  then  the  induced  velocity  at  P,  due  to  the 
disturbance  caused  by  the  SOSL  is  i73 

.  _w*£iiy£i) 


If  P  =  P2,  the  induced  velocity  is 


<X> 
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*(xp2,yp2) 

itPv^l-o2 


it 


~2sin*1 


t]2-o2 

l-o2 


(47) 


Referring  to  Figure  2-1  IB,  the  additional  induced  velocity  inside  the  area  bounded  by  the  tip 
and  the  Mach  line  emanating  from  the  tip  (P  =  P3)  is: 
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Vlx P3>yP3 ) 

Ttp/l^rp 


cos'1 


[o|+Tl2 


(48) 


Again  if  P  =  P4,  the  point  is  out  of  the  zone  of  influence  of  the  SOSL  and  thus  the  induced 
velocity  is  zero. 

The  induced  velocity  at  a  given  point  on  any  wing  geometry  can  now  be  computed  by 
the  proper  superposition  of  the  triangular  SOSL  shown  in  Figures  2-10  and  2-11.  This  is 
because  of  the  linear  nature  of  the  governing  flow-field  Equation  (1).  As  an  example  of 
how  the  above  superposition  principle  works,  consider  the  wing  shown  in  Figure  2-12.  For 
simplicity,  the  slopes  xi  and  arc  constant.  The  wing  AHJD  can  be  represented  by  the 
superposition  of  five  SOSL.  The  first  has  the  planform  AEH  and  source  intensity: 


w(xp,yp)  =  VjCj 


(49) 


where  Xi  is  the  slope  of  the  segment  AB.  The  second  has  the  planform  BIF  and  intensity 


w{xp,yp)  =  (x2~Xi)V- 


(50) 
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MACH  LINE  MACH'UNE 


FIGURE  2-11  A.  TRIANGULAR  SURFACE  SYMMETRIC  ABOUT  X-AXIS 

FOR  SUPERSONIC  SOSL 


FIGURE  2-1  IB.  WING  TIP  EFFECTS  FOR  SUPERSONIC  SOSL 
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and  the  third  the  planform  DJG  and  intensity 

v(xplyp)  =  -%2Vm  (51) 

The  other  two  SOSL  represent  the  tip  effects.  They  are  the  planforms  HJL  and  UL  and  have 
source  intensities  of  opposite  signs  than  those  representing  the  wing. 

The  above  procedure  can  be  applied  to  a  wing  of  general  planform.  The  only 
difference  is  that  for  each  point  in  question,  the  slope  is  not  constant  as  was  the  case  in  the 
simplified  example.  Then  for  some  general  point  located  on  the  wing  surface,  the  total 
induced  velocity  due  to  all  sources  and  sinks  is  found  by  applying  one  of  the  Equations  (43) 
through  (48)  for  each  SOSL.  The  particular  equation  applied  depends  upon  the  location  of 
the  point  relative  to  the  SOSL  and  the  Mach  line  as  discussed  earlier.  These  individual 
contributions  are  then  summed  to  get  the  total  induced  velocity.  Knowing  the  total  induced 
velocity  at  a  point  allows  one  to  calculate  the  pressure  coefficient  at  the  given  point  by 
Equation  (42). 

The  pressure  coefficient  can  be  calculated  at  a  given  number  of  span  wise  and 
chordwise  locations.  The  drag  of  a  given  airfoil  section  at  the  spanwise  station  y  =  yA  is 
then 


Cd  *  ‘clyj  Jo  ^  Cp [X' Ya)  w[x>y*)dx  {52) 

The  total  drag  for  one  fin  of  semispan  b/2  is  then: 

i  fb f* 

Cv=-±[  Cjdy)  dy  (53) 

where  Sw  =  b/2(cr  +  cj.  For  cruciform  fins,  the  total  drag  coefficient  is: 

cDs4-[b/2  Cdc(y)dy  (54) 

swj  o 

If  it  is  desired  to  base  the  drag  coefficient  on  the  body  cross-sectional  area,  the  Equation 
(54)  must  be  multiplied  by  the  factor  S^S^. 
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Equations  (52)  and  (54)  can  be  integrated  by  numerical  quadrature  if  the  generators 
of  the  wing  surface  are  supersonic.  If  the  generators  are  subsonic,  linear  theory  indicates  the 
pressure  coefficients  go  to  infinity  at  the  wing  generators.  Physically,  this  cannot  be  true 
which  means  that  for  a  subsonic  SOSL,  linear  theory  is  not  valid  at  the  SOSL.  The  reason  is 
that  the  velocity  perturbations  in  the  vicinity  of  the  discontinuities  are  no  longer  small, 
violating  one  of  the  assumptions  in  linear  theory.  However,  the  velocity  perturbations  are 
small  a  slight  distance  from  the  SOSL  so  that  linear  theory  can  be  applied.  Numerical 
experiments  indicated  a  distance  of  five  thousandths  of  the  chord  length  from  the  SOSL  is 
sufficient  and  the  value  of  pressure  calculated  at  this  point  can  be  assumed  to  exist  up  to  the 
SOSL. 


The  analysis  using  TDTWT  has  been  illustrated  for  the  axial  force  computation  using 
the  first  term  of  the  boundary  condition  of  Equation  (39).  A  very  similar  process  is  used 
for  the  lift,  roll  and  pitch  damping  computations.  The  reader  is  referred  to  references  74 
through  86  for  the  theoretical  derivations  and  to  Moore  et  al 2,4  for  the  practical  application 
of  the  theories  for  these  force  or  moment  components.  Time  will  not  permit  the  many 
applications  of  TDTWT. 


2. 10  SLENDER  BODY  AND  LINEAR  THEORY  FOR  INTERFERENCE  LIFT 
COMPUTATION*7 


The  method  almost  universally  used  for  including  interference  between  the  various 
missile  components  into  approximate  aeroprediction  codes  is  that  due  to  Pitts,  et  al.*7  There 
are  three  primary  types  of  interference  lift  (note  that  lift  and  normal  force  are  used 
interchangeably  here)  to  be  concerned  with.  These  are  the  effects  on  the  wing  due  to  the 
presence  of  the  body,  the  effect  on  the  body  due  to  the  presence  of  a  wing,  and  finally,  the 
effect  on  an  aft  lifting  surface  due  to  wing  or  body  shed  vortices.  Wing  to  wing  or  shock 
wave  interference  will  not  be  discussed  at  present. 

To  better  understand  the  interference  lift  components,  it  is  instructive  to  examine  the 
total  normal  force  of  a  configuration  as  defined  by  Pitts  et  al.*7  This  is  given  by 

(55) 

+[(*T(fi)  +*S(fl)«]  (  Ctf.)  r+^ATr(v) 

The  first  term  in  Equation  55  is  the  normal  force  of  the  body  alone  including  the 
linear  and  nonlinear  components;  the  second  term  is  the  contribution  of  the  wing  (or  canard) 
including  interference  effects  and  control  deflection;  the  third  term  is  the  contribution  of  the 
tail  including  interference  effects  and  control  deflection;  and  the  last  term  is  the  negative 
downwash  effect  on  the  tail  or  body  due  to  wing  shed  or  body  shed  vortices.  The  K’s 
represent  the  interference  of  the  configuration  with  respect  to  angle  of  attack,  and  the  k’s 
represent  the  interference  with  respect  to  control  deflection.  Each  of  these  interference 
factors  is  estimated2  by  slender  body  or  linear  theory.*7  As  such,  they  are  independent  of 
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angle  of  attack. 

The  various  interference  factors,  as  defined  by  slender  body  theory  (SBT),  are*7: 


kh(b)  ~  2 /it 


(l+r4/s4)[-itan'1-|  (s/i-z/s)  +tc/4 
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Figure  2-13  plots  the  interference  lift  factors  given  by  Equations  (56)  through  (59)  as  a 
function  of  the  body  radius  to  wing  semispan  plus  body  radius  ratio  (r/s). 
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FIGURE  2-13.  SLENDER  BODY  THEORY  INTERFERENCE  LIFT  FACTORS 
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As  the  Mach  number  increases  supersonically,  SBT  gives  values  of 
which  are  too  high  if  the  wing  is  near  the  missile  rear.  This  is  because  much  of  the 
carryover  lift  onto  the  body  is  actually  lost  to  the  wake  of  the  vehicle.  Figure  2-14  illustrates 
this  for  the  no  afterbody,  infinite  afterbody,  and  short  afterbody  cases.  Linear  theory 
formulations  are  available  for  the  infinite  and  no  afterbody  cases  to  replace  Equation  57  if 
the  parameter 

pAJ?{l+X)  [1/  (mP)  +1]  >4  (60) 

Moore2  then  linearly  interpolated  between  the  infinite  and  no  afterbody  cases  as  a  function  of 
the  area  covered  by  the  Mach  lines  to  obtain  for  the  short  afterbody  case. 

Strictly  speaking,  the  methodology  discussed  here  is  limited  to  slender  bodies  with 
triangular  planforms  of  low  aspect  ratio.  Experience  has  shown,  that  if  the  correct  value  of 
wing-alone  lift  is  computed,  the  interference  factors  can  give  very  reasonable  results  for 
wings  which  do  not  have  triangular  planforms  or  even  have  low  aspect  ratio.  Moore2 
showed  how  an  engineering  estimate  of  interference  lift  could  be  obtained,  even  for 
planforms  such  as  that  shown  in  Figure  2-15A.  The  actual  SBT  configuration  is  that  shown 
in  Figure  2-15B.  Since  most  of  the  interference  lift  occurs  near  the  wing  body  juncture, 
reference  (2)  used  approximations  given  by  Equation  (61) 


*  1  +  ^  1  j~l )  G 

(3)^11  =  1  +  (  1^3(3)  J  & 

1  1 1  =  (  1  I  ~  [■ktfCB)  1  j)  & 


(61) 


to  estimate  the  interference  factors  of  the  wing  in  Figure  2- 15 A.  G  in  Equation  (61)  is  the 
ratio  of  the  root  chord  of  the  wing  for  which  the  interference  factor  is  desired  to  that  of  the 
wing  that  slender  body  theory  assumes.  That  is 
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NO  AFTERBODY  (  XAFJ«  0) 


FIGURE  2-14.  DETERMINATION  OF  FOR  HIGH- ASPECT-RATIO 

RANGE  AT  SUPERSONIC  SPEEDS 
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FIGURE  2-15A.  WING  FOR  WHICH  INTERFERENCE  LIFT  IS  DESIRED 


FIGURE  2-15B.  ASSUMED  SLENDER  BODY  REPRESENTATION 
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The  last  two  terms  of  Equation  (55)  are  also  interference  terms.  is  the  lift  on 
the  tail  caused  by  the  vortices  shed  by  the  wing  or  canard  upstream.  CN  is  the  negative 
lift  on  the  afterbody  due  to  wing  shed  vortices.  These  terms  are  also  calculated  analytically 
and  are  given  by: 

_  v(Cl T  [Ku(B)  sin  a  +  kw{B)  sin  6  v]  i  ( sT-iT)  Aw 


WT(V 1 

2 it  (AR)  T(fv-ru)  Aref 

.2  2  2 

Cv  =  -4r 

£w  zw  £  t  rT 

AwVm 

i-  ‘  /4  *4 

Here  i  is  the  tail  interference  factor  given  by  Pitts  et  al87  and  T  is  the  strength  of  the  wing 
shed  vortex. 


2.11  EMPIRICAL  METHODS2  4  6 

It  is  fair  to  wonder  why  approximate  aeroprediction  codes  are  defined  as 
semiempirical  with  all  the  theoretical  methods  discussed  so  far.  The  truth  is,  that  while  these 
methods  allow  the  individual  component  forces  and  moments  to  be  calculated  fairly 
rigorously  at  a  given  Mach  number  or  angle  of  attack,  there  are  still  many  conditions  where 
the  analytical  methods  presented  previously  are  either  not  applicable  or  the  difficulty  in 
applying  then  is  not  worth  the  effort.  In  those  cases,  empirical  methods  are  generally  used. 
The  combination  of  theoretical  and  empirical  techniques  in  a  code  is  thus  why  they  are  called 
semiempirical  codes.  A  few  examples  where  empirical  methods  are  used  are  transonic 
aerodynamics,  body  alone  subsonic  aerodynamics,  rotating  band  or  protuberance 
aerodynamics,  and  base  drag  of  the  body  and  lifting  surfaces.  There  are  actually  analytical 
methods  available  for  transonic  aerodynamic  computations.  However,  most  of  the  methods 
are  inconsistent  from  a  computational  standpoint  with  the  approximate  codes.  What  is  done 
in  many  cases,  is  to  use  the  sophisticated  analytical  tools2,  *• 6  to  estimate  the  transonic 
aerodynamics,  as  a  function  of  key  geometric  parameters,  then  to  include  these  into  an 
engineering  code  in  a  table  lookup  fashion.  Obviously,  for  a  vehicle  that  spends  a  large 
portion  of  its  time  in  the  transonic  flow  region,  0.8  <  M*  <  1.2,  it  would  be  justifiable  to 
use  a  more  sophisticated  estimation  process. 

The  base  drag  empirical  method  will  be  discussed  in  more  detail  in  the  next  section  of 
the  report,  which  deals  with  some  of  the  newer  nonlinear  methods  developed  in  the  past  three 
years. 
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3.0  NEW  APPROXIMATE  AERODYNAMIC  METHODS 


This  part  of  the  report  will  deal  with  many  of  the  new  aerodynamic  prediction 
methods  developed  over  the  past  3  years.  These  methods  include  extension  of  the  SOSET  to 
include  real  gas  effects  (including  two  new  nonlinear  angle-of-attack  pressure  predictors),  an 
improved  version  of  the  Modified  Newtonian  Theory  (IMNT),  and  improvements  to  the 
Allen  and  Perkins  viscous  crossflow  theory;  also  included  are  a  new  nonlinear  wing-alone 
method,  new  nonlinear  wing  body  and  body  wing  interference  methods  due  to  angle  of 
attack,  a  new  nonlinear  wing  body  interference  method  due  to  control  deflection,  a  method 
for  treating  nonlinear  wing  tail  interference,  and  an  improved  base  drag  prediction  model. 

These  new  methods  and  improvements  were  directed  at  three  weak  areas  in  the  NSWCDD 
Aeroprediction  Code  of  1981  (AP81):  (1)  limited  Mach  number  and  inability  to  compute 
temperatures  at  the  surface  for  aeroheadng  calculations,  (2)  lack  of  nonlinear  lift  capability 
except  for  the  body  alone,  and  (3)  base  drag  methodology  that  was  not  robust  enough  in 
terms  of  including  fin  effects. 

3.1  SOSET  EXTENDED  TO  REAL  GASES*9, 90 

The  main  reason  the  fourth  version7  of  the  aeroprediction  code  was  limited  to  Mach 
number  8  was  that,  above  M„=  6  real  gas  effects  start  becoming  important  but,  can  still  be 
neglected  at  M„  =  8.  However,  as  Mach  number  increases  substantially  above  M„  =  6, 
the  need  to  include  real  gas  effects  into  the  aeroprediction  code  increases  if  one  is  interested 
in  inviscid  surface  temperatures.  If  one  is  only  interested  in  forces  and  moments,  real  gas 
effects  have  a  slight  effect  on  the  pitching  moment,  but  only  second-order  effects  on  axial 
and  normal  force®8.  However,  one  of  the  key  issues  in  high-speed  vehicles  is  aerodynamic 
heating,  material  selection,  and  insulation.  Any  excess  weight  can  have  a  strongly  adverse 
impact  on  vehicle  performance.  Thus,  a  simple  yet  accurate  method  of  estimating  vehicle 
surface  temperature  (inviscid)  for  use  in  heat  transfer  analysis  is  needed. 

Figure  3-1 89  is  an  illustration  of  the  importance  of  real  gas  effects.  It  plots  the  static 
temperature  behind  a  normal  shock  for  both  perfect  and  real  gases  at  an  altitude  of  170,000 
ft.  At  this  altitude,  the  speed  of  sound  is  approximately  1100  ft/sec  and  the  freestream  air 
temperature  is  approximately  283  °K.  The  normal  shock  would  occur  in  the  vicinity 
immediately  ahead  of  the  blunted  portion  of  a  seeker  or  the  missile  nose.  Note  that  the 
temperatures  of  interest  to  tactical  weapons  aerodynamicists  can  be  very  high,  for  high  Mach 
number  conditions  assuming  a  perfect  gas.  Also  shown  on  the  figure  are  the  real  gas 
results92.  Note,  in  particular,  the  plot  of  TR/TP,  the  ratio  of  the  real  gas  to  perfect  gas 
temperature.  For  Mach  numbers  of  6  or  less,  this  ratio  is  unity  or  near  unity.  This  is  the 
reason  that  aerodynamic  computations  below  M»  =  6  could  neglect  real  gas  effects  with 
little  error.  However,  as  M.  goes  above  M„  =  6,  the  error  in  temperature  using  the  perfect 
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gas  assumption  becomes  increasingly  large.  This  is  of  particular  importance  to  materials  and 
structures  engineers  designing  the  system  to  withstand  these  temperatures.  Also  shown  in 
Figure  3-1  is  the  melting  point  of  typical  structural  materials  used  in  present-day  missile 
design.  The  actual-use  temperature  is  less  than  the  melting-point  temperature.  For  missiles 
that  fly  at  any  appreciable  time  above  the  maximum-use  temperature  of  a  given  material, 
some  form  of  active  cooling  or  insulation  would  be  required.  This  means  additional  dead 
weight  and,  hence,  less  performance  for  the  missile.  It  is  therefore  obvious  that  a  reasonably 
accurate  estimate  of  temperature  is  essential  for  the  design  of  the  seeker  and  the  structure  of 
the  weapon.  To  meet  the  need  for  a  fairly  accurate  method  of  predicting  surface 
temperature,  SOSET  was  extended  to  include  real  gas  effects.  In  so  doing,  new  approximate 
methods  were  developed  for  angle  of  attack  pressure  prediction  and  an  improved  version  of 
MNT  was  derived.  These  new  methods  will  be  briefly  described. 


STATIC 

TEMPERATURE 
BEHIND 
NORMAL 
SHOCK  T2 
(KxlO-3) 


FIGURE  3-1.  TEMPERATURE  BEHIND  A  NORMAL  SHOCK  AS  A  FUNCTION 
OF  FREES TREAM  MACH  NUMBER  (H  =  170kft) 

SOSET  and  MNT  for  perfect  gases  were  discussed  in  2.1  and  2.3,  respectively. 

Refer  to  2. 1  for  the  SOSET  methodology  and  to  Moore,  et  al.®9, 90  for  the  extension  to  real 
gases.  It  is  noted  that  to  extend  SOSET  to  real  gases  requires  several  things:  (1)  a  cone 
solution  for  real  gases  (pe);  (2)  a  Prandtl-Meyer  Expansion  (PME)  for  real  gases  (pj);  (3)  a 
derivation  of  a  new  pressure  derivative  (dp/ds)2,  where  the  perfect-gas  assumption  has  not 
been  made;  and  (4)  a  way  to  compute  temperature  given  values  of  pressure.*9  After  the  real- 
gas  pressure  derivative  (dp/ds)2  was  derived  and  checked,  it  was  found  that  (dp/ds)2  became 
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negative  for  many  cases,  causing  one  to  choose  between  the  Generalized  Shock  Expansion 
Theory  (GSET  where  jj  =  0)  and  the  tangent  cone  theory  (t?  =  »).  in  comparisons  of  the 
pressure  prediction  to  full  Euler  computations,  it  was  found  that  a  better  way  to  implement 
the  shock  expansion  theory  for  M  ^  6  was  to  redefine  Equation  (10)  as 

P  =  Pe  “  (Pc  -  P2)  hi  (64> 

with  17,  being  an  input  parameter  chosen  by  the  user.  It  was  found  that  a  value  of  ij,  =  0 
gave  slightly  better  pressure  predictions  for  slightly  blunt  configurations,  whereas  a  value  of 
tj,  =  1  gave  better  accuracy  where  bluntness  was  large.  Thus,  final  implementation  of 
SOSET  in  AP93  is  Equation  (64),  with  i*,  as  an  input,  pc  the  rea  *as  tangent  cone  pressure, 
and  P2  the  real-gas  value  of  pressure  computed  from  a  Prandtl-Meyer  expansion. 

To  compute  inviscid  temperatures  (and  other  p'operties)  along  the  surface  of  a  pointed 
or  blunt  body,  the  constancy  of  entropy  along  the  surface  for  perfect,  frozen,  or  equilibrium 
chemically  reacting  flows  is  used.  Knowing  the  value  of  entropy  and  pressure  from  the 
pointed  cone  solution92  or  the  normal  shock  solution  for  a  blunt  body93,  one  can  then  use  the 
thermofit  equations  of  Tannehill  and  Mugge94  and  Srinivasen,  et  al.,95  to  determine  other 
properties,  i.e., 


T  =  Tip,  S) 
P  “  P(P/S) 
a  =  a  ip,  S) 
e  =  eip,S) 


(65) 


The  remaining  properties  at  the  body  surface  can  be  found  from  standard  thermodynamic 
relationships,  i.e., 


h  =  e  +  p/p 


"•  ■  (**) 


r0_  =  constant 


v  =  y/Ti^-h] 

M=V/a 


=  Jif£ 


z=-£— 

pRT 


(66) 
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In  the  process  of  computing  surface  properties,  three  new  pressure  prediction  methods 
were  derived.  The  first  of  these  was  to  give  an  improved  pressure  coefficient  prediction  on 
the  blunt  nose  of  a  missile  configuration  over  that  provided  by  the  MNT.  If  the  pressure 
coefficient  of  MNT  is  defined  as 

(Cp)  Mm  =  £p0sin26e(?  (67) 

then  the  nose  pressure  on  the  blunt  nose  part  of  a  missile  is  given  by 

Cp  =  (Cp)^  -ACp  (68) 


ACP  of  equation  (68)  is  defined  by 

ACp  =  kcosm  (6eg)  [cosfieg-  cos  (6  eg)  J 

where  (6  eq)m  =  25.95  deg,  m  =  2.78,  and 


k  = 


2.416  C„ 


4.606 


2 

0.1507C  ♦ 

Po 


1/2 


1.124 

Ml 


(69) 


Figure  3-2  shows  the  results  of  the  Improved  Modified  Newtonian  theory  (IMNT)  of 
Equations  (68)  and  (69),  compared  to  Equation  (67)  alone,  and  a  full  numerical  solution  of 
the  Euler  equations67  for  a  hemispherical  forebody  at  M.  =  10.  The  IMNT  gives  up  to  7 
percent  improvement  in  pressure  compared  to  the  MNT.  Even  past  the  match  point 
(5eq  <  25.95  deg),  the  IMNT  gives  good  agreement  with  the  numerical  solution  down  to 
6eq  values  of  10  deg.  This  level  of  accuracy  in  pressure  prediction  will  also  translate  into 
more  accurate  drag  computations,  particular  on  bodies  with  large  bluntness. 

The  other  two  pressure  prediction  formulas  have  to  do  with  calculating  the  pressure 
on  a  point  behind  the  blunt  nose  portion  of  the  body  but  at  an  angle  of  attack.  These  are 

Cp(a,<t>)  =  CPmmB  -  (2a)  sin  (20)  cos  (<|>)  +  (Fcos20)a2  + 

(4/3  sin  (20)  cos  (<|>) )  a3  (70) 


where 


(2  -  r 


(1  -  tan20c)  -  (2 


+  ~ )  sin2<|> 


and 


cp  (a,4>) 


'p«. 


(2a )  sin  (20)  cos  (<f>) 
3 


(71) 
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FIGURE  3-2.  SURFACE  PRESSURE  DISTRIBUTION  OVER  A 
HEMISPHERICAL  FOREBODY  AT  M.  =  10 


Equation  (70)  is  used  for  pointed  body  configurations,  as  well  as  for  blunt  body 
configurations  in  the  windward  plane  area  (60°  <<t>  <;  180°).  Equation  (71)  is  used  in  the 
leeward  plane  (<f>  <  60°)  for  configurations  with  blunt  noses.  In  Equation  (70),  (Cp)a  .  0  *s 
the  pressure  coefficient  at  a  =  0,  which  comes  from  Equation  (64).  Figure  3-3  is  an 
example  of  the  application  of  Equation  (70)  to  a  cone  along  with  the  associated  inviscid 
surface  temperatures.  The  approximate  results  are  close  to  the  exact  cone  solution 97 
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FIGURE  3-3.  PERFECT-GAS  COMPARISON  OF  EXACT  AND  APPROXIMATE 

CONE  SOLUTIONS 


Figure  3-4  presents  the  comparison  of  the  present  methodology  for  predicting  in  viscid 
surface  temperatures  on  a  20-percent  blunt  cone  at  a  =  10  deg  and  M„  =  IS.  These 
results  are  compared  to  a  full  numerical  solution  of  the  Euler  equations  (ZEUS)22  for  both 
perfect  and  real  gases.  The  real-gas  temperatures  are  substantially  lower  than  the  perfect-gas 
results  and  also  agree  with  the  full  Euler  solution  except  in  the  vicinity  of  the  overexpansion 
region  past  the  blunt  tip.  Figure  3-3  uses  most  of  the  theory  developed  for  the  approximate 
methodology  in  Equations  (64)  through  (71),  along  with  the  assumptions  used  in  computing 
temperature. 
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AXIAL  LOCATION,  X  (CAL) 


FIGURE  3-4.  COMPARISON  OF  APPROXIMATE  AND  EXACT  TEMPERATURE 
IN  WINDWARD  PLANE  OF  A  20  PERCENT  BLUNT  CONE 
(M  .  =  15,  a  =  10  DEG) 
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3.2  AEROHEATING9* 


The  AP93  methodology  computes  boundary  layer  heating  information  in  the  form  of  a 
heat  transfer  rate,  a  heat  transfer  coefficient,  H;  and  a  recovery  temperature  (adiabatic 
wall  temperature),  Tw,  at  each  computational  point.9®  These  variables  are  related  as  shown 
in  Equation  (72). 


H  = 


(72) 


Tw  is  the  wall  temperature.  For  high-temperature  flows,  the  heat  transfer  coefficient  is  often 
expressed  in  terms  of  enthalpies. 


Hx  = 


(73) 


At  temperatures  above  about  1500°R,  Equation  (73)  is  the  more  rigorously  correct  of 
the  two.  The  heat  transfer  is  normalized  as  shown  in  Equations  (72)  and  (73)  because  the 
coefficients  H  and  H,  remain  fairly  constant  over  a  wide  range  of  wall  temperatures,  even 
though  the  actual  heat  transfer  rate,  may  vary  significantly.  Thus,  since  T„  and  h„  are 
not  functions  of  wall  temperature,  once  a  heating  computation  is  performed  for  a  given  Mach 
number/altitude  combination,  it  need  not  be  repeated  simply  because  of  changes  in  wall 
conditions.  This  weak  coupling  greatly  simplifies  the  problem  of  tracking  the  time-dependent 
thermal  response  of  a  surface  exposed  to  boundary  layer  heating.  The  aerodynamic  solution 
may  be  obtained  first  with  a  code  such  as  AP93,  and  the  results  stored  in  tabular  form  as 
functions  of  a  Mach  number,  altitude,  and  angle  of  attack.  This  information  can  then  be 
accessed  by  an  independent  algorithm  to  compute  the  time-varying  heat  transfer  rates  and  the 
resulting  integrated  surface  temperature  history  along  any  given  trajectory  that  lies  within  the 
limits  of  the  data  matrix. 

The  only  departure  from  the  use  of  true  inviscid  surface  conditions  as  boundary  layer 
edge  properties  occurs  in  the  case  of  blunt  bodies.  The  curvature  of  the  detached  bow 
shocks  associated  with  these  configurations  creates  an  entropy  layer  near  the  body  surface. 
The  inviscid  solution  would  give  a  uniform  boundary  layer  edge  entropy  over  the  entire  body 
equal  to  that  behind  a  normal  shock  at  the  free-stream  Mach  number,  since  this  is  the 
entropy  along  the  inviscid  streamline  that  wets  the  body  surface.  In  reality,  because  of  the 
finite  thickness  of  the  boundary  layer,  the  true  edge  entropy  is  that  which  exists  at  some 
point  in  the  entropy  layer  located  at  a  distance  above  the  surface  equal  to  the  local  boundary 
layer  thickness.  This  entropy  value  is  determined  by  an  iterative  mass  balance  technique.9® 
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Once  appropriate  boundary  layer  edge  conditions  are  determined,  a  series  of 
specialized  analytical  relations  are  used  to  determine  the  aerodynamic  heating  at  various 
locations.  At  the  nose  tip  stagnation  point,  a  simplified  version  of  the  Fay-Riddell  formula" 
gives 


=  0.763Pr-°-6/p^  J 


dv. 

dx 


JV  A~) 


(74) 


The  stagnation  point  velocity  gradient,  dVJdx,  is  determined  from  the  Newtonian 
theory,  assuming  a  spherical  nose  tip.  At  the  nose  tip,  the  flow  will  always  be  laminar. 


If  control  surfaces  are  present,  the  viscous  heating  along  their  leading  edge  stagnation 
lines  is  determined  by  the  Beckwith  and  Gallagher  swept-cylinder  relations100  modified  to 
include  real-gas  effects.101  For  the  laminar  case, 


=  0.57  Pz 


-0.6 


/PoPo 


dV 

~^{haw-hw)  (cosA)1-1 


(75) 


where  A  is  the  leading  edge  sweep  angle  and  dV^dx  is  the  stagnation  line  velocity  gradient 
derived  from  Newtonian  theory,  assuming,  a  cylindrical  leading  edge.  For  turbulent  flow, 


<5r„  t  -  i  .  04Pr-0"6 


(pV) 

(Ho) 


0.8 

o” 


(l^sinA)  °-6 


ter 


(76) 


where  Vp  is  the  flow  velocity  parallel  to  the  leading  edge  stagnation  line  and  the  (*) 
superscript  denotes  evaluation  at  a  reference  enthalpy  given  by102 

A*=0.5  +0.22  (haw-he)  (77) 


The  (e)  subscript  denotes  evaluation  at  the  boundary  layer  edge.  The  laminar  or 
turbulent  status  of  the  flow  is  determined  by  comparison  of  the  Reynolds  number,  based  on 
the  leading  edge  diameter,  to  user-specified  upper  and  lower  limits.  If  Ren  is  below  the 
lower  limit,  laminar  values  are  used.  If  Reu  is  above  the  upper  limit,  fully  turbulent  flow  is 
assumed.  For  intermediate  values  of  Re,, ,  a  linear  combination  of  laminar  and  turbulent 
values  is  computed. 
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For  points  on  the  body,  the  Eckert  reference  enthalpy  flat  plate  formulation  is  used.103 
For  laminar  flow, 


<^=0.332  (Pr*) 


-0.667 


Re* 


(78) 


and  for  the  turbulent  case, 


<2r„(t  =  0.185  (Pr*)'°-6677 


P*V. 


2.564 

In  Re' 

[  ”t\ 

(79) 


N,  and  Nt  are  transformation  factors  that  allow  for  the  approximation  of  three- 
dimensional  (3-D)  effects.  They  are  equal  to  three  and  two,  respectively.  The  laminar  or 
turbulent  flow  character,  is  determined  as  before  by  comparing  the  local  Reynolds  number, 
based  on  boundary  layer  running  length,  to  user-specified  upper  and  lower  limits. 

Heating  rates  on  the  surfaces  of  wings,  fins,  or  canards  are  determined  by  using 
Equations  (78)  and  (79)  but  in  this  case,  N,  and  Nt  are  both  equal  to  one  because  of  the  two- 
dimensional  (2-D)  nature  of  the  flow.  The  degree  of  turbulence  is  determined  in  the  same 
manner  as  for  the  body. 


An  example  of  the  new  aeroheadng  method  is  given  in  Figure  3-5.  Figure  3-5  shows 
the  heat  transfer  rate  on  a  15  degree  half  angle  cone  with  a  nose  radius  of  1.1  inches  as  a 
function  of  distance  along  the  axis  of  symmetry.  Conditions  considered  are  Ma  =  10.6  and 
angle  of  attack  10  degrees.  Comparisons  are  made  with  a  more  complicated  approximate 
technique104  that  uses  streamline  tracking  combined  with  the  axisy metric  analog  to  model  3-D 
effects.  Experimental  data  are  also  shown105  along  with  the  results  from  the  MINIVER101 
code  used  in  a  tangent  cone  mode.  AP  93  and  MINIVER  tend  to  under  predict  the  data  by 
about  10  -  15  percent,  a  performance  that  is  credible  considering  the  simplified  nature  of  the 
solution.  Note  that  the  AP  93  gives  improved  results  over  MINIVER  in  the  vicinity  of  the 
stagnation  region  due  to  the  more  accurate  calculation  of  entropy  at  the  edge  of  the  boundary 
layer  and  more  accurate  real  gas  properties. 


3.3  BASE  DRAG 

The  AP81  estimated  base  drag  using  a  composite  of  empirical  data  for  the  body  alone. 
Also,  an  approximation  was  made  for  the  effect  of  angle-of-attack,  fin  location,  and  fin 
thickness  effects  as  a  function  of  Mach  number  based  on  a  limited  amount  of  data.  As  a 
result,  a  requvSt  was  made  U,  ihc  National  Aeronautics  and  Space  Administration,  Langley 
Research  Center  (NASA/LRC)  to  perform  additional  wind  tunnel  tests,  where  additional  base 
pressure  measurements  could  be  taken  to  try  and  quantify  the  effects  mentioned,  plus  those 
due  to  control  deflection. 
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FIGURE  3-5.  WINDWARD  PLANE  HEAT  TRANSFER  RATES  FOR  1.1-IN. 

NOSE  RADIUS,  15-DEG  HALF-ANGLE  CONE  AT  a  =  10  DEG 

Wilcox  was  the  chief  engineer  for  the  tests  that  were  conducted  and  reported.106, 107 
Eighty-nine  base  pressure  taps  were  placed  around  a  7.2  caliber,  5-inch  diameter  body  with  a 
side  mounted  sting.  These  taps  were  placed  every  22.5  deg  in  circumferential  location  and  at 
several  radii  from  the  body  centroid  toward  the  outer  edge.  The  configuration  matrix  of  data 
taken  is  shown  in  Table  3-1.  The  base  pressure  measured  at  each  of  the  89  orifice  locations 
was  then  averaged  over  its  incremental  base  area  to  get  the  average  base  pressure  at  each 
condition,  of  Table  3-1.  Based  on  these  average  base  pressure  measurements  at  each  test 
condition,  changes  in  base  pressure,  and  hence,  base  drag  because  of  a  particular  physical 
model  change,  or  flight  condition  change  could  be  readily  computed  by  simply  subtracting 
the  two  data  points. 
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TABLE  3-1.  CONFIGURATION  INDEX 


Config 

Fins 

Off 

t/c 

0.05 

0.10 

0.15 

0 

x/c 

1.0 

2.0 

0 

s 

10 

20 

o 

(M.  =  2.0) 

a 

(M.  i  2.5) 

1 

X 

Sweep 

Sweep 

2 

X 

X 

X 

0.5,10 

0 

3 

X 

X 

X 

0.5.10 

0 

4 

X 

X 

X 

0.5,10 

0 

5 

X 

X 

X 

0.5.10 

0 

6 

X 

X 

X 

0,5.10 

0 

7 

X 

X 

X 

0.5.10 

0 

S 

X 

X 

X 

0,5,10 

0 

9 

X 

X 

X 

0.5.10 

0 

10 

X 

X 

X 

0.5.10 

0 

11 

X 

X 

X 

0.5.10 

0 

12 

X 

X 

X 

0.5.10 

0 

13 

X 

X 

X 

0.5.10 

0 

14 

X 

X 

X 

0,5.10 

0 

15 

X 

X 

X 

0.5.10 

0 

16 

X 

X 

X 

0.5.10 

No  data 

Using  the  process  described,  along  with  a  wind  tunnel  data  base  not  available  when 
AP81  was  developed,108  a  new  empirical  estimate  of  base  pressure  coefficient  CPg  was 
derived.  This  new  estimate  is  shown  in  Figure  3-6  and  compared  to  the  AP81  value  of  CPg. 
The  two  curves  are  similar,  with  the  AP93  slightly  higher  than  AP81  for  M.  ^  1.5  and 
slightly  lower  than  AP81  for  M.  ^  3.0.  Body-alone  angle-of-attack  effects  on  base  pressure 
are  then  estimated  by 

(Cp^.-(Cp^..J1+0*011!,J  (80) 

Here,  ( Cp)  WF  a=0  comes  from  Figure  3-6  and  F,,  the  increase  due  to  angle  of  attack  from 
Figure  3-7.  Boattail  and  power-on  effects  on  base  drag  are  estimated  as  present  in  AP81. 
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FIGURE  3-6.  MEAN  BODY-ALONE  BASE  PRESSURE  COEFFICIENT  USED  IN 

AP81  AND  AP93 

At  this  point,  it  is  worth  noting  that,  while  the  databases  of  Moore,  et  al.,  and  Butler, 
et  al.,  helped  to  improve  the  estimate  of  base  pressure  as  a  function  of  Mach  number  and 
angle  of  attack  for  the  body  alone,106- l07- 108  additional  data  are  still  needed  for  a  <  15  deg  at 
all  Mach  numbers.  This  need  is  indicated  by  the  dotted  lines  in  Figures  3-7,  which  are 
extrapolations  from  data  available  for  a  ^  15  deg  and  engineering  judgement.  This  same 
statement  will  also  be  even  more  true  for  fin  effects  due  to  control  deflection  and  angle  of 
attack,  as  will  be  discussed  in  the  following  paragraphs. 

The  total  body  base  pressure  coefficient  for  fins  located  flush  with  the  base  is 


t/c,x/c*  0 


(81) 


where  ( Cpj)  ^,.0 ,  F2,  and  F3  come  from  the  AP93  curve  of  Figures  3-6,  3-8,  and  3-9 
respectively. 
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WILCOX  DATA  IHtF  1081  tXTHAPOUkTION 

OM>2.0  _ _ _ 

0>M*U 
T  M  ■  1.0 
•  M>)J 


A.  (M„  ;>  2) 


B.  (M„  <  2) 


FIGURE  3.7.  PERCENT  INCREASE  IN  BODY-ALONE  BASE  PRESSURE 
COEFFICIENT  DUE  TO  ANGLE  OF  ATTACK 
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FIGURE  3-8.  PERCENT  INCREASE  IN  BASE  PRESSURE  COEFFICIENT  DUE  TO 
COMBINED  EFFECTS  OF  ANGLE  OF  ATTACK  AND  CONTROL  DEFLECTION 

(t/c  »  0) 


DATA  SOURCE  WEf  10B) 


FIGURE  3-9.  ADDITIONAL  PERCENT  INCREASE  IN  BASE  PRESSURE 
COEFFICIENT  DUE  TO  FIN  THICKNESS  AT  VARIOUS  VALUES  OF  |a  +  5| 
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In  Figure  3-8,  no  data  were  taken  for  M*  <  2,106' 107  and  none  could  be  found  in  the 
literature.  Hence,  the  data  for  Moe=  2  are  assumed  to  apply  for  M.  <  2  as  well.  While 
this  is  a  big  assumption,  it  is  believed  to  be  better  than  neglecting  the  base  pressure  effect 
due  to  control  deflection  and  angle  of  attack,  which  other  engineering  aerodynamics  codes 
do.  It  is  also  worth  noting  that  Figure  3-9  indicates  what  is  intuitively  obvious:  for  small 
control  deflections  and  angles  of  attack,  fin  thickness  effects  are  important  in  base  pressure 
estimation,  whereas  for  large  values  of  a  and  6,  the  additional  change  in  CPt  due  to  fin 
thickness  is  minimal. 

The  final  parameter  to  define  the  effect  on  base  pressure  is  fin  location  relative  to  the 
body  base.  This  is  done  through  Equation  (82),  where 


t/c,x/c 


(c’.)„..  *  °-01(Ac^ 


«,6,  t/c,x/c 


(82) 


Here  (  CPb)  ^  m  is  the  body-alone  base  pressure  coefficient  at  a  given  angle  of  attack  given 
by  Equation  (80)  and  ( t/c.x/c  is  the  total  change  due  to  the  presence  of  fins  at 
a  given  a,  5,  t/c,  and  x/c.  An  example  of  ( ACPJ  m  #  tic, x/c  is  given  in  Figure  3-10  for 
M*-  2.0  and  |a  +  5|  =10  deg.  Moore,  et  al.,  showed  other  curves  for  this  parameter.106 
Figure  3-10  shows  that  the  change  in  base  pressure  due  to  all  variables  present  varies  from 
that  at  x/c  =  0,  where  the  fins  dominate  to  that  of  the  body  alone  where  the  fins  have  no 
effect  (x/c  m  2.5). 


FIGURE  3-10.  PERCENT  INCREASE  IN  BASE  PRESSURE  COEFFICIENT  DUE 
TO  FIN  LOCATION  |o  +  5|  =  10  DEG,  M.=  2.0 
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3.4  IMPROVED  METHOD  FOR  BODY-ALONE  NORMAL  FORCE  AND  CENTER 
OF  PRESSURE509  1,0 


The  nom."’ -force  coefficient  of  the  body  alone  is  estimated  by  109 


CN.  +  C| 


(83) 


where  CNl  is  the  linear  term  and  CWj(L  the  nonlinear  term.  The  linear  term  is  predicted  in 
AP81  by  either  SOSET,  second-order  Van  Dyke  combined  with  MNT,  or  empirical 
depending  on  the  Mach  number  range.1,6  The  nonlinear  term  is  estimated  by  the  Allen- 
Perkins  viscous  crossflow  theory.34  No  changes  were  made  in  the  linear  term  of  Equation 
(83)  in  AP93  from  AP81.  Three  changes  in  the  nonlinear  term  of  Equation  (83)  were  made 
for  the  AP93. 


The  nonlinear  term  of  Equation  (83)  is69 


C"«L=TlC<*cSin2a 


(84) 


The  first  change  from  AP81  is  in  the  value  of  ij.  AP81  used  an  incompressible  value 
of  rj  with  no  account  of  compressibility  effects,  although  compressibility  effects  have  been 
clearly  shown.111  The  compressibility  effect  is  shown  in  Figure  3-11 A  along  with  the  line 
drawn  to  represent  the  data.  This  line  is  defined  as 

11  =  (  VtK  +  foz  1.8 

'  8 /  (85) 

q  =  1  fox  MN  >  1.8 


where  q0  is  the  incompressible  value  of  rj  (MN  =  0)  used  in  AP811. 

The  second  change  is  in  the  value  of  the  crossflow  drag  coefficient  used.  This  value 
was  changed  to  allow  the  effect  of  transition  on  the  body  surface  to  affect  the  value  chosen. 
This  affects  the  value  of  for  MN  values  of  0.5  and  less.  Also,  the  value  of  Cdc  is 
slightly  lower  for  0.6  <  MN  <  2.2  than  that  used  in  AP81.  This  is  based  on  the  large 
NASA  Tn-Service  Data  Base.112  The  new  value  of  Cdc  used  in  AP93  is  given  in  Figure  3- 
11B.  If  the  flow  on  the  body  is  a  combination  of  laminar  and  turbulent  (the  case  for  most 
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conditions),  a  value  somewhere  in  between  the  two  values  on  the  Figure  3-1  IB  curve  for 
Mn  <  0.5  will  be  computed.  If  XL  defines  the  length  of  laminar  flow  on  the  body  and  XT 
is  the  total  length,  then  for  MN  $  0.5, 

c“.  =  1-2 -(§j0-8  <86) 


Thus,  if  XL  =  0  so  the  flow  over  the  body  is  fully  turbulent,  a  value  of  Cd^  -  1.2  will  be 
computed,  whereas  a  value  of  0.4  will  be  picked  if  the  flow  is  fully  laminar. 

The  third  change  made  in  AP93  was  in  the  center-of-pressure  location.  AP81  used  a 
weighted  average  of  the  normal  force  center  of  pressure  of  the  linear  term  and  nonlinear 
term,  where  the  nonlinear  term  was  at  the  centroid  of  the  planform  area  in  the 
crossflow  plane  and  the  X^  of  the  linear  term  was  computed  theoretically  or  empirically. 
Both  of  these  values  were  held  constant  as  angle  of  attack  increased,  the  only  change  being 
from  the  changing  values  of  the  normal-force  terms  of  Equation  (83).  In  numerical 
experiments  using  the  NASA  Tri-Service  Missile  Data  Base,  it  was  found  that  the  assumption 
of  a  constant  value  of  center  of  pressure  with  angle  of  attack  was  not  completely  correct.  It 
is  suspected  that  as  angle  of  attack  increases,  the  center  of  pressure  of  the  linear  term  of 
Equation  (83)  changes  and  can  no  longer  be  assumed  to  be  constant.  An  empirical  way  to 
represent  this  change  with  Mach  number  is  given  in  Figure  3-1 1C.  This  change  is  effective 
for  a  ^  10  deg.  Between  a  =  0  and  10  deg,  the  correction  is  implemented  in  a  linear 
fashion  between  zero  at  a  =  0  to  its  full  value  at  a  =  10  deg. 

Figure  3-12  is  an  example  of  the  normal-force  and  center-of-pressure  comparisons  of 
the  AP81,  AP93,  and  experimental  data.  The  data  are  for  a  12.33-caliber  tangent-ogive 
cylinder  configuration  with  a  3. 0-caliber  nose.112  The  improvements  made  in  AP93  give 
significantly  better  results  on  both  CN  and  X^  as  a  function  of  angle  of  attack. 


3.5  WING-ALONE  NONLINEAR  NORMAL  FORCE  AND  CENTER  OF  PRESSURE 

One  of  the  major  reasons  the  AP81  gave  poor  results  at  a  >  10  deg  for  many  missile 
configurations  was  the  failure  to  include  nonlinearities  in  wing  lift.  Using  NASA  and  ONR 
Data  Bases113  114  a  semiempirical  method  was  developed  for  the  nonlinear  wing-alone  normal- 
force  term  analogous  to  the  body-alone  Equations  (83)  and  (84). 109,110  The  nonlinear  term  of 
wing-alone  lift,  therefore,  can  be  defined  as 


f{MN,AR,  A) 


(87) 
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(Mn  =  McoSIN  a) 

FIGURE  3-1 1A.  COMPRESSIBILITY  EFFECTS  ON  CROSSFLOW 
DRAG  PROPORTIONALITY  FACTOR 


FIGURE  3-1  IB.  CROSSFLOW  DRAG  COEFFICL-  -TT  FOR  AN 
OGIVE-CYLINDER  CONFIGURATION 


FIGURE  3-1 1C.  CENTER-OF-PRESSURE  SHIFT  IN  BODY- ALONE 
NORMAL  FORCE  FOR  a  £  10  DEG 
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0  5  10  15  20  25  30 


Angle  of  Attack  (deg) 

FIGURE  3-12.  BODY- ALONE  NORMAL-FORCE  COEFFICIENT  AND 
CENTER  OF  PRESSURE  (  M.  =  3.5) 
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Here,  f(MN,  AR,  X)  is  analogous  to  the  t \Cde  of  the  body  alone  in  Equation  (84).  Since  the 
total  wing-alone  normal  force  is  known  for  a  given  AR,  M.,  X,  and  a,113, 114  and  the  linear 
value  of  lift  is  known  from  the  3-D  thin-wing  theory  or  lifting  surface  theory  from  AP81;  the 
nonlinear  normal  force  of  the  wing  alone  is 


CNml  (Mn,AR,\)  =  Cw  (Mtl,AR,  X)  -  CSl  {Mn,AR,X) 


(88) 


Using  the  data  of  References  113  and  114,  Equation  (88)  values  were  generated  and  a 
parameter  k,  defined  as 


CNaL{Ms,AR,X) 
si  n2a 


(89) 


was  generated.  Tables  of  k,  for  both  high  and  low  Mach  numbers  are  given  in  Tables  3-2 
and  3-3.  The  total  wing-alone  normal  force  in  AP93  is  therefore 

C*w  =  C„l  +  /qsin’a-^-  (90) 

The  second  term  of  Equation  (90)  was  neglected  in  AP81. 

The  center  of  pressure  of  the  wing-alone  lift  was  assumed  to  vary  quadratically 
between  its  linear  theory  value  at  a  =  0  to  the  centroid  of  the  planform  area  (adjusted  for 
thickness  effects)  at  a  =  60  deg. 

Defining  the  center  of  pressure  of  the  wing-alone  linear  term  as  A  and  the  center  of 
pressure  of  the  nonlinear  term  as  B  (both  in  percent  of  mean  geometric  chord),  then  the 
center  of  pressure  of  the  wing  lift  is 


<■ v » * 


A  -  •—  ja  J  {B  -  A]  ♦ 


5400  • 


a*  [A  -  B) 


(91) 


aw  is  the  total  angle  of  attack  in  degrees  on  the  wing.  Figure  3-13  gives  an  example  of  the 
AP93  methodology  compared  to  AP81  and  experimental  data.  This  particular  case  shows 
significant  improvement  in  the  wing-alone  normal  force  of  the  AP93  versus  AP81  when 
compared  to  the  experiment.  However,  no  improvement  in  center  of  pressure  is  obtained 
because  X  =  0  and  the  centroid  of  the  planform  area  is  the  same  as  experimental  data 
suggest. 
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TABLE  3-2.  VALUES  OF  k,  FOR  LOW  MACH  NUMBER 

AR  <  0.5;  Mr  <  4.0 


AR  =  1.0;  M,  <  3.5 


AR  >  2.0;  M*  <  3.5 


TABLE  3-3.  VALUES  OF  k,  FOR  HIGH  MACH  NUMBER 


AR  <  0.5;  M,  >  4.0 


V 

M®sin« 

0.0 

1 - 

0.5 

1.0 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 

a 

5.0 

5.5 

6.0 

CO 

-1.60 

-0.98 

0.23 

0.55 

0.71 

0.82 

0.89 

0.92 

0.95 

0.95 

0.95 

0.95 

0.95 

0.5 

-0.87 

-0.24 

0.33 

0.82 

0.89 

0.92 

0.95 

0.95 

0.95 

0.95 

0.95 

1.0 

-0.31 

0.09 

0.46 

0.68 

0.78 

0.87 

0.91 

0.93 

0.95 

0.95 

0.95 

0.95 

0  95 

AR  =  1.0;  Mr  >  3.5 


XJ 

Moosin® 

0.0 

0.5 

1.0 

D 

2.0 

2.5 

3.0 

3.5 

4.0 

D 

5.0 

5.5 

6.0 

0.0 

-0.39 

-0.39 

-0.29 

0.06 

0.29 

0.48 

0.60 

0.69 

0.75 

0.81 

0.86 

0.91 

0.94 

0.5 

0.14 

0.17 

0.29 

0.46 

0.63 

0.76 

0.85 

0.90 

0.93 

0.95 

0.95 

0.95 

0.95 

1.0 

0.30 

0.50 

0.86 

0.93 

0.94 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 

AR  >  2.0;  Mx  ^  3.5 


V 

Musin'* 

0.0 

0.5 

1.0 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 

H 

5.0 

5.5 

6.0 

0.0 

-0.25 

-0.05 

0.20 

0.50 

0.80 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 

0.5 

0.02 

0.29 

0.80 

0.98 

0.98 

0.97 

0.97 

0.96 

0.95 

0.95 

0.95 

0.95 

0.95 

1.0 

0.66 

1.02 

1.15 

1.18 

1.15 

1  09 

1.02 

0.96 

0.95 

0  95 

0.95 

0.95 

0.95 
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FIGURE  3-13.  WiNG- ALONE  NORMAL-FORCE  COEFFICIENT  AND  CENTER 
OF  PRESSURE  (AR  =  0.5,  X  =  0.0,  M.  =  1.6) 
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3.6  WING-BODY  AND  BODY-WING  NONLINEAR  INTERFERENCE  FACTORS 
DUE  TO  ANGLE  OF  ATTACK  109  110 


The  total  configuration  normal-force  coefficient  at  a  given  angle  of  attack,  control 
deflection  and  Mach  number  is  given  by  Equation  (SS)  repeated  here  for  convenience: 


~  *  [(-^(B)  +  ( *0 ) ®  +  (•ktf(B)  +  ^B(tf))^ 

+[(^r(B)  +  KB(T))a  *  (^T(B)  +  (Cn,)t  *  CnT(^ 


(55) 


Moore,  et  al.,  found  that  the  wing-body  interference  factor  KW(B)  had  the  qualitative 
behavior  as  shown  in  Figure  3-14. 109  At  low  angles  of  attack,  slender-body  theory  appeared 
to  be  a  good  estimate  of  KW(B).  This  estimate  was  adjusted  slightly  for  M.  ^  1.5  by  an 
amount  AKW(B).  At  some  angle  of  attack  defined  as  ae,  KW(B)  seemed  to  decrease  in  a  nearly 
linear  fashion.  The  rate  of  this  decrease  was  a  function  of  Mach  number:  the  higher  the 
Mach  number,  the  larger  the  rate  of  decrease.  At  some  point  defined  as  aD,  the  KW(B) 
appeared  to  reach  a  minimum  and  remain  about  constant.  As  a  result  of  this  analysis,  a 
mathematical  model  was  derived  to  define  KW(B)  in  terms  of  its  slender-body  theory  value 
[KW(b)1sb  and  an  empirical  correction  derived  from  several  databases. u2- n3- 1,4  This  model 
given  in  Figure  3-14  is 

*«»  ’  [*.(»]»  *  t0r 


KK(B)  ~  [*V(B)]« 


■  [AK«»].-o  *  -^r(#  -  ec)  (f3)  foz  (92; 

»  i 
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FIGURE  3-14.  QUALITATIVE  BEHAVIOR  OF  WING-BODY  INTERFERENCE 
FACTORS  AS  A  FUNCTION  OF  ANGLE  OF  ATTACK 
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The  empirical  corrections  to  KW(B)  are  also  in  a  form  that  can  be  defined 
mathematically  as  opposed  to  a  table  lookup  procedure.  These  equations  for 


[A-^W<B> 


dKt 


W{B) 


da 


a  r,  a 


are  as  follows: 


P*™.]...'0-22 


for  Af.sl.o 


[A^(B)]a.o  =  -°-44[w--1-5] 


[A**<B>U  =  ° 


for  1.0<Mm£l.5 
for  Mm>  1.5 


=  -(0.00283W.  +  0.025) 

da 


&S2.Q 


a  c  =  12.5  -  1 . 0SMm  -2.5  for  AffsO .  5 

a  c  =  12.5  -  6 . 25Af„  for  AR=1 . 0 


«c=4.5  <■  2.25M.  -  2.25Mf“£ff™|j’j;  for  AS*2 . 0 

Mm>2.0 
a  c=0 


a  D  =  33.3  -  8. 19AT.  +  0 . 82Mi 
a  D  =  25.3  -  6.62  Mm  *  0.66  Ml 


for  A.  =  0 

for  k  =  1 . 0 
for  0  <  k < 1 . 0 


(93) 


(94) 


(95) 
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The  semi-empirical  model  for  was  also  defined  in  terms  of  its  slender  body  or 
linear  theory  value,  plus  a  correction  due  to  nonlinearities  associated  with  angle  of  attack. 
The  mathematical  model  for  was  defined  as  109 


K, 


•fl(» o 


=  + 


r/s 

0.5 


[A*B(V)]«.q 


+ 


^1*1  • 


(97) 


Unfortunately,  a  mathematical  model  for  [AKB(W)]a.0  and  dHC^l/da  was  difficult  to  define 
because  of  the  variability  of  the  constants  as  a  function  of  the  parameters  of  interest.  As  a 
result,  a  three-parameter  table  lookup  for  these  two  parameters  is  used  in  AP93  based  on  the 
data  in  Table  3-4.  The  parameters  in  the  table  lookup  include  M»,  X,  and  AR.  Linear 
interpolation  is  used. 

In  Equations  (92)  and  (96),  the  factor 

r/s 

0.5 


appears.  This  is  because  the  NASA  Tri-Service  Missile  Data  Base  is  based  on  r/s  =  0.5, 
and  Pitts,  et  al.,  indicates  that  the  aerodynamics  vary  linearly  with  r/s.*7  This  assumption  is 
inherent  in  the  semiempirical  models  for  KW(B)  and  K^. 
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TABLE  3-4.  DATA  FOR  BODY-WING  NONLINEAR  SEMIEMPIRIC AL 

INTERFERENCE  MODEL 


Data  for  [ 

«-0 

Mach  Number 

Aspect 

Ratio 

Taper 

Ratio 

£  0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

4.5 

£  0.25 

0.  0.5, 
1.0 

-0.1 

-0.1 

0.5 

0.6 

0.7 

0.8 

0.7 

0.5 

0.3 

0.5 

0.5 

-0.28 

-0.1 

0.13 

0.11 

0.05 

-0.02 

-0.06 

0 

0 

1.0 

0.5 

-0.26 

-0.2 

0.15 

0.21 

0.15 

0 

0 

0 

0 

^  2.0 

0.5 

-0.13 

-0.04 

0.12 

0.43 

-0.16 

0 

0.37 

-0.08 

-0.16 

0.5 

0 

-0.3 

0.06 

0.26 

0.28 

0.17 

0.12 

0.14 

0 

0 

2  2.0 

0 

-0.2 

0.1 

0.12 

0.52 

0.12 

0.15 

0.22 

-0.06 

-0.22 

0.5 

1.0 

-0.16 

0.08 

0.26 

0.14 

-0.12 

0 

-0.05 

-0.10 

0 

St  2.0 

1.0 

-0.2 

0.1 

0.12 

0.45 

-0.02 

0.11 

0.28 

-0.17 

-0.3 

Data  for  dfK^  ]/da 

Mach  Number 

Aspect 

Ratio 

Taper 

Ratio 

£  0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

4.5 

S!  0.25 

0,  0.5 
1.0 

0.018 

0.013 

-0.010 

-0.023 

-0.013 

-0.022 

•0.031 

-0.025 

-0.031 

0.5 

0.5 

0.019 

0.010 

-0.008 

•0.010 

-0.013 

-0.013 

-0.013 

-0.012 

-0.012 

1.0 

0.5 

0.013 

0.010 

-0.007 

-0.013 

-0.020 

•0.017 

-0.012 

-0.012 

-0.012 

2:  2.0 

0.5 

0.010 

0.011 

0 

-0.013 

-0.010 

-0.017 

-0.040 

-0.012 

-0.012 

0.5 

0 

0.033 

0.022 

0 

•0.007 

•0.010 

-0.008 

-0,014 

-0.012 

-0.012 

2:  2.0 

0 

0.010 

0.010 

-0.007 

-0.020 

-0.011 

-0.020 

-0.023 

-0.012 

-0.012 

0.5 

1.0 

0.019 

0 

-0.019 

-0.010 

-0.007 

-0.013 

-0.014 

-0.012 

-0.012 

2:  2.0 

1.0 

0.010 

0.01 

-0.007 

-0.017 

0 

-0.017 

-0.026 

-0.012 

-0.012 
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In  examining  cases  where  r/s  is  small,  it  was  found  that  at  high  angles  of  attack,  the 
wing-alone  solution  was  not  recovered  properly  through  the  process,  Equations  (92)  and  (97) 
To  remedy  this  situation,  the  AP93  nonlinear  interference  factors  were  blended  into  those 
predicted  by  slender-body  or  linear  theory  as  r/s  became  small.  The  specific  equations  used 
to  do  this  are 


For  r/s  >  0.25 

KWB)  =  [K*(B)]AP93 

KBM  =  [^bmLlpsj 


(98a) 


For  0.05  <  r/s  <  0.25 

*«»  1*«»>  W-(  (r/s-°  •  05)  /0 . 2 

(98b) 

kb(w)  =[-KB(io]sflr-/  [-^s ( »o ]ssr “[■Kfl ( io Japj j)  (r/ s_0  •  05)  /0 . 2 

LT  \  LT  J 


For  r/s  <  0.05 

KM(B)  ~  [**<b,]5BT  >  KB{U)  -  Fb(w]bbt  (98c) 

In  essence,  the  model  represented  by  Equations  (98a)  through  (98c)  uses  the  nonlinear 
interference  factors  for  r/s  values  greater  than  0.25;  they  use  a  blend  of  slender-body  or 
linear  theory  and  the  nonlinear  values  of  interference  factors  for  r/s  values  between  0.05  and 
0.25.  They  also  use  the  slender-body  or  linear  theory  values  for  r/s  values  less  than  0.05. 
Hence,  when  the  body  vanishes  (r/s  =  0),  the  wing-alone  solution  will  be  automatically 
recovered  in  a  smoother  and  more  accurate  way. 

Figure  3-15  is  an  example  of  the  normal  force  on  the  wing  in  the  presence  of  the 
body  and  the  normal  force  on  the  body  in  the  presence  of  the  wing  using  AP93  theory,  the 
AP81  theory,  and  compared  to  experimental  data.  Note  that 


CN*tB)  ~  CNWKH(B) 
CNBtm  =  CNmKB  CP0 


(99) 
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Angto  of  Attack  (d«g) 


FIGURE  3-15.  WING-BODY  AND  BODY-WING  INTERFERENCE  AS  A 
FUNCTION  OF  a  (AR  =  2.0,  X  =  0,  M.  =  1.2) 

Hence,  Figure  3-15  is  actually  a  representation  of  the  normal-force  coefficient  on  the  wing 
and  additional  normal  force  on  the  body  due  to  the  wing.  Thus,  Equation  (99)  is  a 
representation  of  the  accuracy  of  not  only  KW(B)  and  K^,  but  CNk  in  conjunction  with  the 
interference  factors.  This  is  a  more  true  indication  of  the  accuracy  of  the  code  because  there 
are  actually  two  of  the  component  force  terms  that  make  up  Equation  (39).  As  seen  in 
Figure  3-15,  the  AP93  methodology  is  clearly  superior  to  the  AP81  theory  as  angle  of  attack 
increases. 

The  center  of  pressure  of  the  new  value  of  normal  force  of  the  wing  in  the  presence 
of  the  body  estimated  by  Equation  (92)  is  assumed  to  remain  at  the  values  of  the  wing-alone 
solution  of  AP93  given  by  Equation  (91).  The  center  of  pressure  of  the  additional  lift  on  the 
body  due  to  the  presence  of  the  wing  is  estimated  using  the  AP81  method,  which  is  either 
slender-body  or  linearized  theory.  These  values  are  modified  for  short  afterbodies.2 


3-30 


NSWCDD/TR-93/55 1 


In  exercising  the  AP93  on  missile  configurations  in  the  transonic  speed  regime 
(0.6  <  M  <  2.0),  it  was  found  that  some  of  the  nonlinear  lift  associated  with  small  aspect 
ratio  fins  (AR  <  1.4)  was  lost  due  to  shock-wave  formation.  An  empirical  approach  in  the 
AP81  accounted  for  a  certain  amount  of  linear  lift  loss.  This  appeared  to  be  satisfactory  for 
the  larger  aspect  ratio  fins,  where  the  nonlinear  normal-force  term  with  angle  of  attack  was 
negative.  However,  when  the  fins  have  a  positive  nonlinear  normal  force  due  to  angle  of 
attack,  some  of  this  force  appears  to  be  lost  with  shock  waves.  This  loss  was  estimated 
empirically  as  a  function  of  Mach  number  and  angle  of  attack  for  a  wing  that  had  an  area- 
to-body  reference  area  of  about  one.  These  data  for  ACN  losses  due  to  compressibility 
effects  are  given  in  Table  3-5.  A  two  parameter  linear  interpolation  is  made  from  Table  3-5 
for  a  given  M.  and  a  to  compute  ACN  .  ACN  is  further  degraded  for  taper  ratio  for  values 
of  X  <  0.5.  The  specific  equations  for  ACN  are 


'ref 


A 

=  -0.2A  C, 


for  A*0.5 


for  0 . 1  £  A  £  0 . 5 


for  A  s  0 . 1 


(100) 


TABLE  3-5.  LOSS  OF  WING  NONLINEAR  NORMAL  FORCE  DUE  TO  SHOCK- 
WAVE  EFFECTS  IN  TRANSONIC  FLOW 


a  +  5 

,  deg| 

M. 

0 

5 

10 

15 

20 

25 

30 

35 

2s  40 

S  0.4 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

O.0000 

0.0000 

0.0000 

0.0000 

MM 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0220 

-0.2060 

-0.6890 

-0.9500 

-1.300 

0.8 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0531 

-0.2200 

-0.7100 

-1.010 

-1.400 

1.2 

0.0000 

0.0000 

-0.0093 

-0.0293 

■0.1651 

-.04167 

-0.7629 

-1.070 

-1.500 

1.5 

0.0000 

0.0000 

-0.0653 

-0.1111 

-0.1556 

-0.4444 

-0.7000 

-1.070 

-1.500 

2.0 

0.0000 

0.0000 

-0.0076 

-0.0376 

-0.1502 

-0.1142 

•0.0951 

-0.0700 

0.0500 

S2.5 

0.0000 

0.0000 

0.0000 

0.000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 
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3.7  NONLINEAR  WING-BODY  INTERFERENCE  FACTOR  DUE  TO  CONTROL 
DEFLECTION* 


Initially,  it  was  planned  to  use  slender-body  theory  for  the  interference  factors  kW(B) 
and  ,  as  currently  done  in  AP81.  This  plan  was  based  on  results  comparing 
computations  (using  Equations  (55)  where  all  the  nonlinearities  are  included)  with 
experimental  data  at  6  =  0  for  both  body-tail  and  wing-body-tail  or  dorsal-body-tail 
configurations.109  These  comparisons  were  good  and  seemed  to  indicate  that  new  technology 
was  superior  to  existing  engineering  approaches.  However,  when  results  were  examined  for 
configurations  that  had  control  deflections  on  either  the  aft  or  forward  lifting  surface,  they 
were  found  to  be  not  as  good  as  desired.  This  led  to  the  conclusion  that  nonlinear 
interference  factors,  due  to  control  deflection,  were  also  required  to  improve  the  performance 
of  AP93  when  compared  to  experimental  data. 

The  approach  taken  was  to  use  the  AP93  with  the  non-linearities  of  wing-alone,  wing- 
body,  and  body-wing  interference  effects  due  to  angle  of  attack  included,  use  the  slender- 
body  estimates  of  kW(B)  and  k^  for  control  deflection,  and  derive  empirical  modifications  to 
kW(B)  based  on  numerical  experiments  compared  to  actual  missile  data.  Because  kW(B) 
appears  in  the  vortex  lift  on  the  tail  due  to  canard  or  wing  shed  vortices,  the  numerical 
experiments  were  conducted  with  canard  body-tail  configurations. 

Referring  to  Equation  (55),  the  vortex  normal-force  coefficient  on  the  tail  is*7 


'N, 


TiV) 


t  [^w<b)  sina  +  FkWiB)  sin 5 ^i(sT  ~  rr)Av 
2  71  ( AR )  jfff,  -  TfjjAr0t 


(101) 


Equation  (101)  has  a  factor  F  that  multiplies  the  term  due  to  control  deflection  in  the 
wing-tail  vortex  lift.  This  factor  is  needed  in  addition  to  the  nonlinearity  for  kW(B),  partly 
because  the  negative  afterbody  lift  due  to  control  deflection  is  not  presently  modeled  in  either 
AP81  or  AP93.  This  term  is  defined  by  Equation  (63). 

The  main  reason  this  term  was  not  included  in  the  AP81  code  was  that  it  required  an 
estimate  of  fT,  which  is  the  position  of  the  canard  shed  vortex  at  the  tail.  Also,  Nielsen,  et 
al.,  indicated  that  this  term  was  generally  much  smaller  than  that  computed  by  Equation 
(101). 114  To  account  for  this  term,  a  vortex  tracking  algorithm  or  an  empirical  correction  to 
the  term  in  Equation  (101)  is  needed.  For  angles  of  attack  much  greater  than  25  or  30  deg, 
a  vortex  tracking  algorithm  may  be  needed.  However,  up  to  a  of  about  30  deg,  a  nonlinear 
model  of  interference  effects  resulting  from  control  deflection  was  developed  by  defining 
kW(B)  as  a  function  of  angle  of  attack  and  Mach  number  and  F  as  a  function  of  Mach  number 
and  angle  of  attack. 


Using  the  work  of  Nielsen,  et  al.,  McKinney,  and  Smith,  et  al.,  for  low  Mach 
number,114,  n5, 1,7  a  semiempirical  nonlinear  model  for  kW(B)  and  the  parameter  F  were 
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derived  from  numerical  experiments.  The  mathematical  model  for  kw(B)  is  based  on  slender- 
body  theory  similar  to  kW(B)  and  k^  and  modified  for  angle  of  attack  or  control  deflection. 
In  general,  it  was  found  that 

* W{B )  =  *"1  +  <-2(lai«’l  ' 

F  =  C3(W,  |a„|)  (102) 

More  specifically,  kW(B),  C,,  Cj,  and  F  are  defmed  in  Figure  3-16  for  Mach  numbers  where 
data  are  available.  For  Mach  numbers  less  than  0.8  and  greater  than  4.6,  the  equations 
derived  for  those  conditions  have  been  used.  The  current  method  for  using  the  empirical 
estimate  for  kW(B)  from  Figure  3-16  is  to  linearly  interpolate  between  Mach  numbers  for  a 
given  value  of  a,  5,  and  M„. 

The  model  in  Figure  3-16  has  a  lot  of  similarities  to  the  nonlinear  KW(B)  model  already 
discussed:  at  low  angle  of  attack,  slender-body  theory  gives  a  reasonable  estimate  of  kw(B). 
However,  as  angle  of  attack  increases,  kW(B)  decreases  up  to  low  supersonic  Mach  numbers. 
For  higher  supersonic  Mach  numbers,  kW(B)  actually  increases  at  higher  angles  of  attack, 
presumably  due  to  compressibility  effects.  Also,  for  low  angles  of  attack,  a  value  of  F  near 
one  is  found  for  the  vortex  lift  model,  indicating  again  reasonable  accuracy  of  the  theory  in 
reference  87.  However,  as  angle  of  attack  is  increased,  F  increases  above  one  for  many 
Mach  numbers.  That  is,  Equation  (101)  gives  values  of  csnV)  too  small  due  to  control 
deflection  of  a  forward  surface.  As  already  mentioned,  this  is  most  probably  due  to  the 
neglect  of  the  effect  on  the  afterbody  Equation  (63),  which  accounts  for  a  greater  percentage 
of  the  afterbody  effect  compared  to  the  Equation  (101)  results,  as  angle  of  attack  increases. 
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F=  1.1 


MS. 8 


lfl«wl  <  24.0  — ky^g)  -  1-4[kw(mlsB 

If  |ocw|  >  24.0  — ►kwfBj  =  1.4  [.000794  l*w|2  -  .0933  1**1  +  2.71] 


M  —  1.1 


If  |0CWI  <  15.0— »kvy(B)  -  1-3[kw(B)lsB 

lfl«wl  >  15.0  —  kw(B)  =  1.3  [.00087  |aw|2  -  .0825  l«wl  + 1.98] 


M  —  1.5 

If  |awl  s  10.0  —  kyy(B)  =  .9[kw(B)]sB 

If  l«wl  >  10.0  —  k^B)  =  .9[kw(B)]SB  -  .01S[I«J  -  10.0] 

If  l«w|  s  20.0  —  F  =  .8 

If  l*wl  >  20.0  —  F=  .8  +  ,10[|ocw|  -  20.0] 

M  =  2.0 

If  l«wl  s  10.0  — ky^B)  =  -9[kW(8)]sB 

If  l«wl  >  10.0  — kv^B,  =  .9[kw(B)]sB  -  .005[|acwl  -  10.0] 

If  l«wl  <  20.0  —  F  =  .8 

If  l«wl  >  20.0  —  F  =  .8  +  -17[|otw|  -  20.0] 

M  =  2.3 

If  |ocw|  S  20.0  — ky^g)  =  -9[kw(B)]sB 

I1"1*1  >  20.0  — kv^B)  =  .91^8^8  -  .005(l*wl  -  20.0] 

If  |ocw|  S  30.0  — F  =  .9 

If  l«wl  >  30.0  —  F  =  .9  +  .15[l<xw|  -  30.0] 

M  =  2.87 

If  l«wl  S  20.0  — ky^B)  =  -9[kw(8)]sB 

\l  l«wl  >  20.0  —  kyy(B)  =  .9[kw(B)lsB  ~  -005[|ac  |  -  20.0] 

If  l«wl  <  30.0  —  F  =  .9 

If  lotyyl  >  30.0  —  F  =  .9  +  ,17[l«w|  -  30.0] 

M  =  3.95 

^w(B)  =  -8[k  w(B)lsB 

If  l«wl  s  40.0 -F  =  0.9 

If  l«wl  >  40.0  —  F  =  0.9  +  ,4[|ocw|  -  40.0] 

M>4,6 

l«wl  s  20.0  — kyy(B)  =  0.75[kw(B)]sB 
f  |0CW|  >  20.0  —  kyy(B)  =  0.75[kyy(B)]sB  +  .01[l«wl  ~  20.0] 
If  l«wl  s  35.0  — F  =  .9 
If  l<xw|  >  35.0  — F  =  .9  +  ,3[l<xw|  -  35.0] 

where  *w  =  *  +  6 


FIGURE  3-16.  NONLINEAR  WING-BODY  INTERFERENCE  MODEL  DUE 

TO  CONTROL  DEFLECTION 
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4.0  SUMMARY  OF  METHODS  IN  1993  VERSION  OF  NSWCDD  AEROPREDICTION 
CODE  (AP93)  AND  COMPARISON  WITH  EXPERIMENT4,  “■ 

The  methods  used  for  computing  forces  and  moments  in  the  AP93  are  summarized  in 
Tables  4-1,  4-2,  and  4-3.  Note  that  the  code  can  now  be  useful  for  computing  aerothermal 
information  as  well  as  forces  and  moments.  This  means  the  code  now  has  five  uses: 

a.  Providing  inputs  to  flight  dynamics  models  that  estimate  range  or  miss  distance 

b.  Assessing  static  stability  of  various  missile  configurations 

c.  Assessing  various  design  parameters  in  terms  of  optimizing  the  configuration 

d.  Assessing  structural  integrity  using  the  loads  portion  of  the  code 

e.  Assessing  aerothermal  aspects  of  a  design  using  the  heat  transfer  coefficients  at 
high  Mach  numbers. 

As  seen  in  Tables  4-1,  4-2,  and  4-3,  there  are  many  methods  that  go  into  the  overall  makeup  of 
a  component  build  up  code,  such  as  the  APC.  The  past  20  years  have  shown  that  this  type  of 
code  can  be  quite  useful  when  used  in  preliminary  or  conceptual  design  studies  to  provide  down 
selection  on  many  configuration  alternatives  in  a  fairly  accurate  and  cost-effective  manner. 

Most  of  the  methods  listed  in  the  tables  have  been  briefly  summarized  in  sections  of  the  report. 

Several  different  complete  missile  configurations  have  been  considered  in  the  validation 
of  the  AP93  code  compared  to  experimental  data.*- 109  A  sample  of  several  of  the  flight 
conditions  on  a  few  of  the  configurations  considered  will  be  given  here.  Also,  there  will  be 
comparisons  with  AP81  or  other  SOTA  aeroprediction  codes  when  such  results  are  available  in 
the  literature.  Funds  were  not  available  to  do  a  thorough  comparison. 

The  first  case  for  comparison  of  the  AP93  and  AP81  is  the  configuration  shown  in 
Figure  4-1  A.  The  body  shown  has  a  three-caliber  tangent  ogive  nose  with  a  total  length  of 
12.33  calibers.  It  has  aspect  ratio  2.0  tails  and  0.1  dorsals.  Mach  numbers  of  4.5  and  10  are 
considered,  and  comparisons  are  made  with  the  ZEUS  code.  Results  of  these  comparisons  in 
terms  of  normal  force  coefficient  and  center  of  pressure  as  a  function  of  angle  of  attack  are 
shown  in  Figure  4- IB.  Center  of  pressure  results  show  the  AP93  within  two  percent  of  the 
body  length  compared  to  the  ZEUS  computations  at  all  angles  of  attack  considered.  On  the 
other  hand,  the  AP81  center  of  pressure  results  differ  by  as  much  as  8  percent  of  the  body 
length  from  the  ZEUS  code.  In  examining  normal  force  coefficient  comparisons,  it  is  seen  that 
at  Mach  4.5  the  AP93  is  within  5  percent  of  the  ZEUS  code,  whereas  the  AP81  results  arc  low 
by  as  much  as  30  percent  due  to  the  omission  of  nonlinear  wing-alone  and  interference  lift.  At 
M  =  10,  the  normal  force  of  AP93  is  within  13  percent  of  the  ZEUS  code,  whereas  the  AP81 
results  are  off  by  as  much  as  40  percent. 
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TABLE  4-1.  AP93  METHODS  FOR  BODY-ALONE  AERODYNAMICS 


Component/ 

Mach  Nunber  Rerion 

Subsonic 
M„  <  OR 

■  |  1  ■ 

Low  Supersonic 
I.:<M„S2.4 

High  Supersonic 

2.4  <  M«,  <  6  0 

Hypersonic 
>  6.0 

Nose  Wave  Drag 

- 

Scmicmpirical  based 
on  Ruler  Solutions 

iiiii 

SOSET  plus  IMNT 

SOSET  plus  IMNT 
Modified  for  Real  Gases 

- 

Wu  and  Aoyoma 

Second-Order  Van  Dyke 

SOSET 

SOSET  for  Real  Cases 

Skin  Friction  Drat; 

Van  Driest  II 

Base  Drac 

Improveil  Empirical  Method 

Acrohcaling  Information 

- 

- 

- 

- 

SOSET  plus  IMNT  for 
Real  Oases 

Inviscid  Lift  ami 

Pitching  Moment 

Empirical 

Scmicmpirical  based 
on  Ruler  Solutions 

Tsien  First-Order 
Crossflow 

SOSET 

SOSET  for  Real  Cases 

Viscous  Lift  and  Pitch 
Moment 

Improved  Allen  and  Pet  kins  Crossflow 
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TABLE  4-2.  AP93  METHODS  FOR  WING-ALONE  AND  INTERFERENCE 

AERODYNAMICS 


Component/ 
Mach  Number  Rcrinn 


Subsonic 
M„  <  O.S 


Transonic 

O  S  <  <  1.2 

1  \  MSB 

High  Supersonic 

2  .4  <  M„  S  6.0 

Hypersonic 

M_  >  fi.O 

Empirical 

Linear  Theory  plus 
MNT 

n  (SE) 

S 

SE  plus  MNT  for  Real 
Gases  Along  Strips 

Van  Driest  11 

Empirical 


Improved  Empirical 


3DTNVT 

Empirical 


Slcndcr-Body  Theory  or  Linear  Theory  Modified  for  Short  Afterbodies 

Empirical 


Slcndcr-Dody  Theory 

Empirical 


Line  Vortex  Theory  with  Empirical  Modifications  fur  k\y(B)Tcrm  and  Nonlincarilics 


None  Present 


SE  plus  MNT  foe  Real 
Gases 
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TABLE  4-3.  AP93  METHODS  FOR  DYNAMIC  DERIVATIVES 


Componcnt/Mach 
Number  Rcsinn 

Subsonic 

M«  <  0.8 

Transonic 

0.8  SMm<  1.2 

Low  Supersonic 
l.2SM«S2.4 

High  Supersonic  Hypersonic 

2.4  <  M« «;  6.0  M„  >  6.0 

Body  Alone 

Empirical 

Wing  and  Interference 
Roll  Damnine  Moment 

Lifting  Surface  Tlicory 

Empirical 

Linear  Thin  Wing  Theory 
or  Strip  Theory 

Wins  Maenus  Moment 

Assumed  Zero 

Wing  and  Interference 
Pitch  Damnine  Moment 

Lifing  Surface  Theory 

Empirical 

I  ES 

Linear  Thin  Wing  Theory 
or  Strip  Theory 
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FIGURE  4-1  A.  BODY-DORSAL-TAIL  CONFIGURATION  USED  FOR 
COMPARING  ZEUS,  IAP,  AND  OAP  COMPUTATIONS 
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FIGURE  4- IB.  COMPARISON  OF  PRESENT  NORMAL  FORCE  COEFFICIENT  AND 
CENTER  OF  PRESSURE  COMPUTATIONS  WITH  THE  ZEUS  CODE  FOR  THE 
DORSAL-BODY-TAIL  CONFIGURATION  OF  FIGURE  4-1 
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The  second  configuration,  Figure  4-2A,  is  taken  from  Howard  and  Dunn.118  The 
dorsals  have  an  aspect  ratio  of  0.12  and  tail  surfaces  have  an  aspect  ratio  of  4.  The 
aeroprediction  code  will  not  handle  the  configuration  as  shown  at  the  top  of  Figure  4-2A. 
Experience  has  shown  it  necessary  to  keep  the  lifting  surface  area,  centroid  of  area,  span, 
taper  ratio,  and  aspect  ratio  the  same  in  the  configuration  modification  process.  This  means  the 
tip  and  root  chord  of  the  dorsal  and  tail  surfaces  had  to  be  adjusted  with  these  constraints  in 
mind.  The  new  adjusted  configuration  is  shown  at  the  bottom  of  Figure  4-2A.  Hence,  this 
configuration  has  all  parameters  outside  the  empirical  data  base  for  use  in  the  AP93  including 
Mach  number,  aspect  ratio,  body  configuration,  and  r  /  s. 

Howard  and  Dunn  showed  only  normal  force  coefficient  results  for  the  body-tail  and 
body-dorsal-tail  configurations  at  M  =  0.1. 1,1  Results  of  the  AP81,  AP93,  and  Missile 
DATCOM  are  shown  in  Figure  4-2B  compared  to  experiment  for  both  the  body-tail  and  body- 
dorsal-tail  configurations.  For  the  wing-body  case,  the  AP93,  and  Missile  DATCOM  produce 
almost  identical  results;  both  show  higher  CN  values  than  experiment,  particularly  at  low  angles 
of  attack.  It  is  not  clear  why  this  discrepancy  exists.  The  AP81  results,  which  have  the  older 
values  of  C*  and  no  nonlinear  wing  lift,  show  even  higher  results  than  either  the  AP93  or 
Missile  DATCOM. 

The  body-dorsal-tail  configuration  results  of  Figure  4-2B  show  that  the  AP93  is  clearly 
superior  to  both  the  AP81  and  Missile  DATCOM.  Normal  force  errors  of  the  AP93  are  less 
than  5  percent  at  all  conditions,  whereas  errors  of  the  AP81  and  Missile  DATCOM  are  as  high 
as  40  and  50  percent,  respectively.  The  fundamental  reason  for  the  AP93  success  is  the 
nonlinear  wing-alone  normal  force  and  interference  factor  methodology.  At  a  =  30°,  the 
body-dorsal  and  dorsal-body  contributes  about  %  of  the  total  configuration  normal  force. 

The  third  configuration  for  validation  of  the  new  semiempirical  methodology  is  shown  in 
Figure  4-3A.  This  configuration  also  differs  substantially  from  the  geometry  characteristics 
from  whic»  die  new  semiempirical  methodology  was  derived.  The  body  is  21.2  versus  12.33 
calibers  long  with  a  2-caliber  Von  Kantian  versus  a  3-caliber  tangent-ogive  nose.  The  dorsals 
and  tail  surfaces  have  aspect  ratios  of  0.36  and  2.14,  respectively,  both  at  the  outer  edge  of  the 
data  base. 

Wind  tunnel  data  exist  for  both  the  body-tail  and  body-dorsal-tail  configuration  for  Mach 
numbers  of  2.3  to  4.6  and  at  several  roll  orientations.119  Comparisons  are  made  at  <t>  =  0° 
roll  and  at  Mach  numbers  of  2.3  and  4.6  for  both  the  body-tail  and  body-dorsal-tail 
configurations.  Results  of  these  comparisons  are  shown  in  Figure  4-3B  for  the  body-tail  and 
Figure  4-3C  for  the  body-dorsal-tail.  The  AP93  results  are  within  the  expected  accuracy 
bounds  on  normal  force,  center  of  pressure,  and  pitching  moment.  While  AP81  results  are  not 
shown  for  clarity,  significant  improvements  in  normal  force  for  both  body-tail  and  body-dorsal- 
tail  configurations  occur  with  less  significant  improvements  in  center  of  pressure.  As  noted  in 
the  comparisons,  the  AP93  is  slightly  superior  to  Missile  3  for  most  pitching  moments22  and  the 
two  codes  (AP93  and  Missile  3)  are  about  equal  in  normal  force  prediction. 
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CONFIGURATION  TESTED  IN  WIND  TUNNEL  (FROM 
REFERENCE  29  WHERE  DIMENSIONS  ARE  IN  INCHES) 


MODIFIED  CONFIGURATION  USED  IN  AEROPREDICTION  COMPUTATIONS 


PARAMETERSFOR  BOTH  MODELS 


(AR)t=4.0 

(AR)d  =  .12 


bt=3.76  in.  XT=16  (ALe)t=24°  Aj=3.54in.2 

bo  =  1.32  in.  AD=77  (ALe)o=60°  AD  =  14.2in.2 


FIGURE  4-2A.  CONFIGURATION  USED  FOR  COMPARISON  WITH 
MISSILE  DATCOM  AND  EXPERIMENT 
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FIGURE  4-2B.  COMPARISON  OF  PRESENT  NORMAL  FORCE  COEFFICIENT 
WITH  THAT  PREDICTED  BY  MISSILE  DATCOM  AND  EXPERIMENT 
FOR  CONFIGURATION  OF  FIGURE  4-2A 


4-9 


NSWCDD/TR-93/55 1 


1.75 


.50  / 

0 1 


k  .07  AT  .50  LOCATION 


.25 


T  I 


1.75 


-H  w-  -25 


SECTION  A- A 


SECTION  B-B 


FIGURE  4-3 A.  DORSAL-BODY-TAIL  CONFIGURATION  USED  FOR 
COMPARING  MISSILE  3,  AP93,  AND  AP81  COMPUTATIONS 
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FIGURE  4-3B.  COMPARISON  OF  PRESENT  NORMAL  FORCE  COEFFICIENT  AND  PITCHING  MOMENT 
COEFFICIENTS  WITH  MISSILE  3  ON  THE  CONFIGURATION  OF  FIGURE  4-3A  (BODY-TAIL  PORTION  OF) 
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FIGURE  4-3C.  COMPARISON  OF  PRESENT  NORMAL  FORCE  COEFFICIENT  AND 
PITCHING  MOMENT  COEFFICIENTS  WITH  MISSILE  3  ON  CONFIGURATION  OF  4-3 
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A  fourth  case  considered  is  the  canard-body-tail  care  shown  in  Figure  4-4 A. 120  The 
configuration  is  somewhat  of  an  extreme  case  for  the  body-alone  aerodynamics  because  it  is  a 
hundred  percent  blunt  and  is  about  22.3  calibers  long.  The  configuration  tested  in  the  wind 
tunnel  has  hangers  attached  to  the  body  for  aircraft  carry  and  launch.  However,  tests  were 
conducted  with  and  without  the  hangers,  and  the  results  showed  that  CN  and  CM  were 
unchanged  but  CA  was  increased  with  the  hangers  present.  The  AP93  and  AP81  theoretical 
computations  are  compared  to  the  corrected  data  of  Groves  and  Fournier,120  where  the  hangers 
have  been  omitted.  Results  are  given  in  Figures  4-4B  through  4-41  for  Mach  numbers  of  0.8, 
2.86,  and  4.63  and  at  canard  deflections  of  0,  10,  and  20  deg.  Examining  Figures  4-4B 
through  4-41,  it  is  shown  that  AP93  gives  good  agreement  with  experimental  data  under  almost 
all  conditions.  Significant  improvements  of  the  AP93  over  the  AP81  are  seen  at  the  lower 
Mach  numbers  and  at  the  higher  Mach  number,  higher  angle-of-attack  conditions. 


FIGURE  4-4 A.  CANARD-BODY-TAIL  CONFIGURATION  USED  IN 

VALIDATION  PROCESS120 

In  analyzing  why  this  improvement  occurs  at  those  conditions,  it  is  noted  that  the  aspect 
ratio  of  the  tail  surfaces  of  the  configuration  of  Figure  4-4A  is  about  0.87  and  that  of  the 
canard  is  about  1.7.  Examining  Tables  4-2  and  4-3,  the  nonlinearity  in  wing-alone  lift  is  small 
for  Mach  numbers  greater  than  about  1.5.  As  normal  Mach  number  increases,  [M.  sin 
(a  +  5)]  and  Mach  numbers  exceed  about  3.5  to  4.0,  nonlinearity  due  to  compressibility 
becomes  important.  As  long  as  the  aerodynamics  are  fairly  linear,  the  AP81  gives  good 
results  up  to  moderate  angles  of  attack.  However,  when  nonlinearities  are  present,  the  AP93 
shows  significant  improvement.  This  improvement  is  the  greatest  on  the  Figure  4-4A 
configuration  at  low  Mach  number  because  the  nonlinear  normal-force  term  on  the  canards  is 
negative,  whereas  that  of  the  tails  is  positive.  The  combination  produces  a  strong  couple  in 
terms  of  the  pitching  moment  as  evidenced  by  Figures  4-4A  through  4-41.  A  good  nonlinear 
capability,  such  as  that  present  in  the  AP93,  is  absolutely  essential  to  get  accurate  stability  and 
control  information  for  these  cases.  Just  examining  Figure  4-4B,  the  center  of  pressure  of  the 
AP81  at  a  =  20  deg  differs  from  the  experimental  data  by  -9.4  percent  of  the  body  length 
versus  1.3  percent  for  the  AP93. 
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Angt*  ol  Attack  (6«g) 


FIGURE  4-4B.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION  OF  FIGURE  4-4A 

(M.  =  2.86,  6  =  0°) 
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Angta  ot  Attack  (dag) 


FIGURE  4-4C.  NORMAL- AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION  OF  FIGURE  4-4A 

(M.  =  4.63,  5  =  0°) 
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FIGURE  4-4D.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION  OF  FIGURE  4-4A 

(M.  =  0.8,  5  =  10°) 
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FIGURE  4-4E.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION  OF  FIGURE  4-4A 

(M.  =  2.86,  5  =  10°) 
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Angt*  04  Attack  (dag) 


FIGURE  4-4F.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION  OF  FIGURE  4-4A 

(M.  =  4.63,  5  =  10°) 
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FIGURE  4-4G.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION  OF  FIGURE  4-4A 

(M.  =  0.8,  5  =  20°) 
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FIGURE  4-4H.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION  OF  FIGURE  4-4A 

(M.  =  2.86,  8  -  20°) 
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FIGURE  4-41.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION  OF  FIGURE  4-4A 

(M.  =  4.63,  5  =  20°) 


4-21 


NSWCDD/TR-93/55 1 


A  fifth  case  considered  in  the  validation  of  the  AP93  code  is  a  configuration 
representative  of  the  SPARROW  missile  tested  at  NASA/LRC.115, 116  The  configuration  tested 
and  reported  by  Monta  is  shown  in  Figure  4-5A.116  The  configuration  tested  by  McKinney  is 
just  like  the  one  tested  by  Monta,  except  it  had  wiring  tunnels  and  wave  guides  present.115 
These  appendages  add  to  the  normal  force  and  pitching  moment,  but  were  not  accounted  for  in 
the  analytical  computations  that  are  presented  in  Figure  4-5.  The  Monta  configuration  did  not 
have  these  appendages  present  and  was  the  main  set  of  data  used  for  the  nonlinear  empirical 
model  validation.  These  results  are  distinguished  in  Figure  4-5  by  the  fact  that  the  cases  that 
had  wave  guides  present  are  indicated. 

Results  of  the  AP81  and  AP93,  compared  to  the  experiment  for  the  configuration  of 
Figure  4-5A,  are  shown  in  Figure  4-5B  through  4-5G.  Results  are  presented  in  terms  of  CN 
and  CM  versus  angle  of  attack  for  various  control  deflections  and  Mach  numbers.  The 
nonlinear  models  with  and  without  control  deflection  show  the  AP93  code  agreeing  much  closer 
to  the  data  at  all  Mach  numbers  than  the  linearized  approaches  of  AP81.  On  the  other  hand, 
the  fact  that  the  body-alone  normal  force  of  AP81  had  the  nonlinearities  included  makes  the 
comparisons  to  experimental  data  better  than  it  would  be  otherwise. 

In  examing  Figure  4-5B,  it  is  seen  that  both  CN  and  CM  of  AP93  agree  with  the 
experiment  at  8=  0  and  8  =  10  deg  for  M*  =  1.5  whereas,  CN  and  CM  of  the  AP81  are 
both  considerably  in  error  as  angle  of  attack  increases  above  5  to  10  deg.  For  M.  =  2.35 
(Figure  4-5C),  both  CN  and  CM  of  AP  93  at  5=  0  and  20  deg  agree  with  the  data.  Again, 
AP81  yields  considerable  error  at  a  ^  10  deg,  although  the  error  is  decreasing  with  increasing 
Mach  number.  For  M„  =  3.95  (Figure  4-5C),  AP81  gives  acceptable  results  for  CN  and  CM 
up  to  a  =  15  to  20  deg  and  at  both  5  =  0  or  20  deg.  The  comparison  with  data  gets  worse 
above  a=  20  deg,  whereas  AP93  comparisons  show  good  agreement  at  all  values  of  a  and  8. 
The  same  statements  basically  hold  true  for  the  M*  =  4.6  comparisons  (Figure  4-5C). 

Figures  4-5F  and  4-5G  show  the  comparisons  of  AP81  and  AP93  to  the  McKinney 
data,115  which  is  the  same  configuration  as  that  of  Figure  4-5 A,  except  that  wave  guides  and 
wiring  tunnels  were  attached  to  the  wind  tunnel  model.  As  already  mentioned,  no  account  was 
taken  for  these  appendages  in  the  analytical  computations.  Note  that  AP93  agrees  much  more 
with  the  data  than  AP81  for  both  M„  =  2.3  and  4.6  at  all  values  of  8.  In  comparing  the  wind 
tunnel  data  for  the  cases  with  and  without  appendages,  it  can  be  seen  that  the  appendages  add 
only  a  few  percent  to  the  aerodynamics. 
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FIGURE  4-5A.  AIR-TO-AIR  MISSILE  CONFIGURATION  USED  IN 
VALIDATION  PROCESS42-43 


2.35,  Delta -0  Mach-2.35,  Delta 


FIGURE  4-5C.  NORMAL-FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  FOR 
CONFIGURATION  OF  FIGURE  4-5A  FOR  VARIOUS  MACH  NUMBERS  AND 
CONTROL  DEFLECTIONS  (CONTINUED) 


Mech-3.95,  Delta *0  Mech-3.95,  Delta 
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FIGURE  4-5D.  NORMAL-FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  FOR 
CONFIGURATION  OF  FIGURE  4-5A  FOR  VARIOUS  MACH  NUMBERS  AND 
CONTROL  DEFLECTIONS  (CONTINUED) 
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FIGURE  4-5E.  NORMAL-FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  FOR 
CONFIGURATION  OF  FIGURE  4-5A  FOR  VARIOUSMACH  NUMBERS 
AND  CONTROL  DEFLECTIONS  (CONTINUED) 


Mach-2.3,  Della—  20  Mach-2.3.  Delta— 12 
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FIGURE  4-5F.  NORMAL-FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  FOR 
CONFIGURATION  OF  FIGURE  4-5A  FOR  VARIOUS  MACH  NUMBERS 
AND  CONTROL  DEFLECTIONS  (CONTINUED) 


4.6,  Data— 20  Mach -4.6.  Data— 12 
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FIGURE  4-5G.  NORMAL-FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  FOR 
CONFIGURATION  OF  FIGURE  4-5A  FOR  VARIOUS  MACH  NUMBERS 
AND  CONTROL  DEFLECTIONS  (CONTINUED) 
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A  sixth  and  final  case  used  in  the  validation  and  development  of  the  nonlinear 
aerodynamics  model  is  shown  in  Figure  4-6A.  Note  that  in  Figure  4-6A,  two  configurations 
were  actually  tested,  one  that  had  a  full-tail  surface  and  a  second  that  had  a  partial  cutout 
removed.117  The  AP93  will  not  handle  the  partial-wing  configuration  as  it  stands,  so  an 
engineering  model  of  this  wing  must  be  created.  Experience  has  shown  that  the  lifting  surface 
area,  aspect  ratio,  span,  leading  edge  sweep  angle,  and  centroid  of  the  presented  area,  must  be 
held  constant.  The  chord  is  varied  so  as  to  meet  these  constraints.  Hence,  the  configuration 
that  represents  the  partial-wing  results  is  the  body  canard  of  Figure  4-6A,  plus  the  AP93 
representation  of  the  partial  tail  shown  in  the  lower  right  of  Figure  4-6A. 
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FIGURE  4-6A.  CANARD-CONTROLLED  MISSILE  CONFIGURATION  WITH 
FULL-TAIL,  PARTIAL-TAIL,  AND  AP93  REPRESENTATION  OF  PARTIAL  TAIL 

FOR  USE  IN  VALIDATION  PROCESS117 


Figures  4-6B  through  4-6D  present  the  comparison  of  the  AP93  with  wind  tunnel  test 
data.  Data  were  only  available  at  M.  =  0.2;  however,  this  complements  the  previous  data  set 
for  the  SPARROW  missile  in  the  sense  that  no  subsonic  data  were  available  for  that  case.  The 
full-tail  and  partial-tail  results  are  denoted  on  the  figure.  Some  results  were  available  from 
reference  117  for  the  Missile  Datcom.17  These  results  are  also  shown  where  available. 
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COMPARISON  OF  AP93  TO  WIND  TUNNEL  DATA  AND  MISSILE 
7OR  NORMAL-FORCE  AND  PITCHING  MOMEN  COEFFICIENTS 
FIGURE  4-6A  CONFIGURATION  (M_  =  0.2,  5C  =  -20°) 
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D.  COMPARISON  OF  AP93  TO  WIND  TUNNEL  DATA  AND  MISSILE 
[  FOR  NORMAL-FORCE  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  FIGURE  4-6A  CONFIGURATION  (M.  =  0.2,  5  =  20°) 
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As  seen  in  the  figure,  the  AP93  gives  improved  results  for  pitching  moment  and  normal  force 
for  most  conditions,  compared  to  the  Missile  Datacom.  While  center  of  pressure  is  not  shown, 
the  AP93  computations  are  generally  within  the  goal  of  ±  4  percent  of  the  body  length.  For 
example,  at  a  =  30  deg,  6  =  -20  deg,  x^  for  the  data,  AP93  and  Missile  Datacom  are  5.39, 
4.91,  and  3.75  calibers,  respectively,  with  respect  to  the  moment  reference  point.  This 
represents  errors  of  2.1  and  7.3  percent  of  the  body  length,  respectively,  for  the  AP93  and 
Missile  Datacom  codes. 

Many  other  cases  have  also  been  considered  in  the  validation  of  the  new  AP93  code.*’109 
In  general,  it  has  been  found  that,  on  average,  the  AP93  code  has  reduced  the  normal  force  and 
center  of  pressure  errors  of  the  AP81  code  by  half,  and  reduced  the  axial  force  errors  by  about 
twenty-five  percent.  There  are  cases  where  AP81  actually  does  better  than  AP93.  However, 
these  are  quite  rare,  and  in  averaging  several  hundred  data  points  for  various  configurations,  at 
various  Mach  numbers  and,  at  5°  increments  in  angle  of  attack  from  0  to  30°,  the  reduction  in 
errors  of  AP93  over  AP81  is  significant.  While  no  equivalent  systematic  comparison  with 
other  SOTA  codes  has  been  made,  the  AP93  was  superior  to  other  engineering  codes  at  most 
conditions  where  comparisons  were  made. 
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5.0  SUMMARY 


In  summary,  various  types  of  aeroprediction  codes  have  been  reviewed.  Some  of  the 
more  important  conventional  approximate  aerodynamic  techniques  have  been  summarized. 
Recent  new  methods  applicable  to  the  semiempirical  class  of  aeroprediction  codes  have  been 
given  in  slightly  more  detail.  Six  complete  missile  configurations  have  been  used  for 
comparing  the  NSWCDD  AP93  code  to  experiment  and  other  approximate  codes.  It  was 
seen  that  this  new  nonlinear  theory  gives  significant  improvement  in  aerodynamic  estimation 
compared  to  the  AP81. 

Some  areas  where  new  technology  is  still  needed  for  semiempirical  codes  is  in  skin 
friction  drag  treatment  with  angle  of  attack,  transonic  aerodynamics,  base  drag  at  a  >15 
deg.,  and  a  simple,  accurate  method  for  estimating  roll-dependent  aerodynamics. 
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8.0  SYMBOLS  AND  DEFINITIONS 


Planform  area  of  the  body  or  wing  in  the  crossflow  plane  (ft2) 

Reference  area  (maximum  cross-sectional  area  of  body  if  a  body  is  present 
or  planform  area  of  wing  if  wing-alone)  (fit2) 

Planform  area  of  wing  in  crossflow  plane  (ft2) 

Speed  of  sound  (ft/sec) 

Aspect  ratio  =  b2/Aw 

Wing  span  (not  including  body)  (ft) 

Total,  base,  and  skin  friction  axial  force  coefficients  respectively 

Drag  Coefficient  -  — DTaJf— 

Vip  mvlAt9f 

Crossflow  drag  coefficient 

Mean  skin  friction  coefficient  based  on  freestream  Reynolds  number  (R*). 

Pitching  moment  coefficient  (based  on  reference  area  and  body  diameter  if 
body  present  or  mean  aerodynamic  chord  if  wing  alone) 

Spanwise  pitching  moment  of  wing  airfoil  section 

Pitch  damping  moment  coefficient  derivative 

Normal  Force  Coefficient  (  NozwalForce ) 

Vap  „vlAzef 

Spanwise  normal  force  of  wing  airfoil  section 
Body  alone  normal  force  coefficient 

Negative  afterbody  normal-force  coefficient  due  to  canard  or  wing  shed 
vortices 

Additional  normal-force  coefficient  on  body  due  to  presence  of  wing 
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AC, 


N, 


am 


'Nr 


'Nm 


'AT, 


r(v) 


' N . 


H(B) 


AC, 


N, 


"(B) 


(Cpb)nf,  „ 
(Cpb)  a.  6, t/c.x/c 


Cr 

Ct 


d 


d«f 

e 


F 


F„  F2,  F3 


Additional  normal-force  coefficient  on  body  due  to  a  control  deflection  of 
the  wing 

Linear  component  of  normal-force  coefficient 
Nonlinear  component  of  normal-force  coefficient 

Negative  normal-force  coefficient  component  on  tail  due  to  wing  or  canard 
shed  vortex 

Normal-force  coefficient  of  wing  in  presence  of  body 

Additional  normal-force  coefficient  of  wing  in  presence  of  body  due  to  a 
wing  deflection 

Normal-force  coefficient  derivative 

.  P~Pm 

Pressure  Coefficient  ( - - ) 

Vfep.vi 

Base  pressure  coefficient 

Base  pressure  coefficient  with  no  fins  present  and  at  angle  of  attack 

Base  pressure  coefficient  with  fins  present  of  some  t/c,  x/c,  8,  and  body  at 
some  a 

Stagnation  pressure  coefficient 
Root  chord  (ft) 

Tip  chord  (ft) 

Body  diameter  (ft) 

Reference  body  diameter  (ft) 

Internal  energy  (ft2/sec2) 

Dimensionless  empirical  factor  used  in  tail  normal-force  coefficient  term 
due  to  canard  or  wing  shed  vortices  to  approximate  nonlinear  effects  due  to 
a  control  deflection 

Symbols  defining  parameters  used  in  base  drag  empirical  model 


fw,  f,  Lateral  location  of  wing  or  tail  vortex  (measured  in  feet  from  body  center 

line) 


8-2 


NSWCDD/TR-93/:51 


H  Heat  transfer  coefficient  based  on  wall  local  temperature  (ft-lb)/(ft2-sec-°R) 

Ho  Total  enthalpy  (ft2/sec2) 

H,  Heat  transfer  coefficient  based  on  wall  local  specific  enthalpy 

[slug/(ft2-sec)] 

h  Specific  enthalpy  (ft2/sec2) 

h^  Adiabatic  wall  specific  enthalpy  (ft2/ sec2) 

1^  Specific  enthalpy  at  outer  edge  of  boundary  layer  (ft2/sec2) 

hT  Height  of  wing  or  canard  shed  vortex  at  tail  center  of  pressure  (ft) 

h^  Specific  enthalpy  at  wall  (ft2/sec2) 

h*  Reference  value  of  specific  enthalpy  (ft2/sec2) 

!  Tail  interference  factor 

Empirical  factor  defined  in  wing-alone  nonlinear  normal-force  coefficient 
term 

KB(W)  Ratio  of  additional  body  normal-force  coefficient  derivative  due  to  presence 

of  wing  to  wing-alone  normal-force  coefficient  derivative  at  5  =  0  deg 

KW(B)  Ratio  of  normal-force  coefficient  derivative  of  wing  in  presence  of  body  to 

that  of  wing  alone  at  5  =  0  deg 

kB(W)  Ratio  of  additional  body  normal-force  coefficient  derivative  due  to  presence 

of  wing  at  a  control  deflection  to  that  of  the  wing  alone  at  a  =  0 

kW(B)  Ratio  of  wing  normal-force  coefficient  derivative  in  presence  of  body  due 

to  a  control  deflection  to  that  of  wing  alone  at  a  -A  0  deg 

[kwreJsB  Value  of  kW(B)  calculated  by  slender-body  theory  at  a  =  0 

AKB(W)  ,  AKW(B)  Nonlinear  corrections  to  KB(W)  and  KW(B)  due  to  angle  of  attack 

1  Length  (ft) 

lN  Nose  length  (can  be  in  calibers  or  feet) 

LT  Linear  Theory 

M  Mach  number  =  V/a 
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Mn  Normal  Mach  number  to  body  axis  =  M  sin  a 

N,,  Nt  Transformation  factors  used  in  Eckert  reference  enthalpy  to  approximate 

three-dimensional  effects  for  laminar  and  turbulent  flow 
(  =  3  and  2,  respectively) 

p  Pressure  (lb/ft2)  or  roll  rate  (rad/ sec) 

pc  Pressure  of  a  cone  of  given  half  angle  (lb/ft2) 

Pr  Prandtl  number 

<%  Pitch  Rate  (rad/sec) 

<£„  Heat  transfer  rate  (ft-lb)/(ft2-sec)  at  wall 

qWi Heat  transfer  rate  at  wall  for  laminar  or  turbulent  flow,  respectively 
R  Gas  constant  [  for  air  R  =  1716  ft-lb/(slug  -  °R)] 

Re  Reynolds  Number  =  -2-^ 

P 

(Re),  Critical  Reynolds  number  where  flow  transitions  from  laminar  to  turbulent 

flow 

ReD  Reynolds  number  based  on  diameter  of  wing  leading  edge  bluntness 

r  Radius  of  body  (ft) 

r„  Radius  of  nose  tip  (ft) 

rw  r,  Radius  of  body  at  wing  or  tail  locations 

r/s  Ratio  of  body  radius  to  wing  or  tail  semispan  plus  the  body  radius 

S  Entropy  (ft-lb)/(slug  -  °Rankine) 

s  Distance  along  body  surface  in  SOSET  (also  wing  or  tail  semispan  plus  the 

body  radius  in  wing-body  lift  methodology) 

SB  Slender -body  theory 

T  Temperature  (°R  or  °K) 

TW,T0,  Tw  Adiabatic  vail,  total,  and  wall  temperature,  respectively 
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Tail  thickness  to  its  root  chord 
Tail  thickness  to  body  diameter 
Perturbation  velocity  components,  (ft/sec) 

Velocity  (ft/sec) 

Velocity  at  edge  of  boundary  layer  (ft/sec) 

Velocity  parallel  to  leading  edge  of  wing  (ft/sec) 

Distance  along  the  axis  of  symmetry  measured  positive  aft  of  nose  tip  (feet 
or  calibers) 

Parameter  used  in  base  drag  methodology  to  represent  the  number  of  chord 
lengths  from  the  base  (measured  positive  upstream  of  base) 

Center  of  pressure  (in  feet  or  calibers  from  some  reference  point  that  can 
be  specified) 

Laminar  and  turbulent  flow  lengths  on  body  (ft) 

S panwise  center  of  pressure  of  wing  semispan 
Compressibility  factor 
Angle  of  attack  (degrees) 

Rate  of  change  of  angle  of  attack  (deg/sec) 

Angle  of  attack  where  wing-body  interference  factor  starts  decreasing  from 
its  slender-body  theory  value  (degrees) 

Angle  of  attack  where  the  wing-body  interference  factor  reaches  a 
minimum  (degrees) 


Local  angle  of  attack  of  wing  or  tail  (  a  +  5W  or  a  +  respectively,  in 
degrees) 
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0  Jm2-1  or  Jl-M2  depending  on  whether  flow  is  supersonic  or  subsonic. 

Also,  Mach  angle,  P=sin'1  (1/M)  . 

b  Control  deflection  (degrees) 

Seq  Angle  between  a  tangent  to  the  body  surface  at  a  given  point  and  the 

velocity  vector  (degrees) 

5W,  <$t  Deflection  of  wing  or  tail  surfaces  (degrees),  positive  leading  edge  up 

0  Velocity  potential 

0  Circumferential  position  around  body  where  0  =  0  is  leeward  plane 

(degrees) 

A  Taper  ratio  of  a  lifting  surface  =  c,/cr 

0,  ,f,  First  order  axial  and  crossflow  solutions  of  velocity  potential  equation 

*2  Second  order  particular  solution  to  full  potential  equation 

t/  Parameter  used  in  SOSET  and  also  used  in  viscous  crossflow  theory  for 

nonlinear  body  normal  force  (in  this  context,  it  is  the  normal  force  of  a 
circular  cylinder  of  given  length-to-diameter  ratio  to  that  of  a  cylinder  of 
infinite  length) 

Vo  Value  of  r}  in  viscous  crossflow  theory  for  MN  =  0 

Mo  <M*  Viscosity  coefficient  at  stagnation  or  reference  conditions,  respectively 

(slug/ft-sec) 

p,  p0  ,  p*  Density  of  air  at  local,  stagnation,  or  reference  conditions,  respectively 

(slugs/ft3) 

y  Specific  heat  ratio 

6  Local  body  slope  at  a  given  point  (degrees) 

0C  Cone  half  angle 

A  Leading  edge  sweep  angle  of  wing  or  tail  (degrees) 

oo  Free-stream  conditions 

2- D  Two  dimensional 

3- D  Three  dimensional 
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3DTWT 

AP81 

AP93 

APC 

BD 

BL 

FNS 

GSET 

IMNT 

MNT 

NASA/LRL 

NS 

NSWCDD 

PNS 

SE 

SOSET 

SOTA 

TAT 


3-D  thin  wing  theory 

Aeroprediction  1981 

Aeroprediction  1993 

Aeroprediction  code 

Base  Drag 

Boundary  Layer 

Full  Navier-Stokes 

Generalized  shock-expansion  theory 

Improved  modified  Newtonian  theory 

Modified  Newtonian  theory 

National  Aeronautics  and  Space  Administration/Langley  Research  Center 
Navier-Stokes 

Naval  Surface  Warfare  Center,  Dahlgren  Division 
Parabolized  Navier-Stokes 
Shock  expansion 

Second-order  shock-expansion  theory 

State-of-the-art 

Turn-Around  Time 


TLNS 


Thin  Layer  Navier-Stokes 
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